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Abstract 

We analyze the chiral symmetry breaking of QCD in the Coulomb gauge. Using 
flow equations, we derive the renormalized gap equation and the Bethe-Salpeter 
equation and show that they are finite in both UV and IR regions. No additional UV 
renormalization is required in the chiral limit. We take into account the hyperfine 
interaction as well as chiral symmetry breaking and obtain the tt — p mass splitting 
caused by the instantaneous and dynamical interactions. 


1 Introduction 

The connection between the fundamental quantum chromodynamics (QCD) and the 
constituent quark model (CQM) is still a puzzling feature in hadron physics. While the 
CQM appears to be based on a rather simple vacuum, the QCD has a complicate non- 
perturbative vacuum structure manifested by the color confinement and the dynamical 
breaking of the chiral symmetry. In this respect, one of the most intriguing phenomena 
in hadron physics may be the mass splitting between the pion (spin 0) and the p meson 
(spin 1), which seems to exibit an interplay between the spin-spin interactions and the 
chiral symmetry breaking. 

In the CQM the spin hyperfine interaction, fj'^Tctxjj'ijj'^Tctxjj, which is a nonrelativistic 
reduction of one gluon exchange is believed to produce the tt — p mass splitting |^. 
However, in the nonrelativistic constituent quark model one obtains the smaller value 
of mass splitting than the experimental observation. Also the typical behavior of the 
hyperfine interaction yields the splitting as as/mgrug and restricts the applications of 
the nonrelativistic CQM to mesons containing at least one heavy quark. Alternatively, 
the 7T — p mass splitting is explained in QCD by the dynamical chiral symmetry breaking 
caused by the chiral noninvariant QCD vacuum. According to the ’t Hooft anomaly 
condition |Q, the chiral symmetry in QCD is broken in the Nambu-Goldstone mode and 
the pion is the Goldstone boson. However, this argument is kinematical and the detailed 
dynamical mechanism of chiral symmetry breaking in QCD remains unclear. 

Even before the existence of QCD, Nambu and Jona-Lasinio 0 successfully described 
the dynamical breaking of chiral symmetry by the fermion pair condensation. In QCD 
a number of approaches have been applied to the problem of chiral symmetry break¬ 
ing. Lattice gauge simulations, though limited to light quarks, predict a nonzero quark 
condensate and show a nonlinear dependence on the light quark mass. In continuum a 
dynamical mass is obtained by solving some sort of Dyson-Schwinger equation in QCD 
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[Q . Moreover the chiral perturbation theory provides the results deviated from the exact 
chiral limit where Goldstone mode is given by = 0. 

Our approach referred as the pairing or the BCS vacuum model is close to the original 
suggestion by Nambu and Jona-Lasinio, which was exploited later by number of authors 

0 ], §• Given a chiral-invariant interaction between massless quarks, '0lT'0'0lT'0, 

one obtains a gap equation by seeking a lowest-energy solution, i.e. the ground vacuum 
state which turns out to be noninvariant under chiral symmetry. Subsequent many-body 
solution of the Bethe-Salpeter equation in this vacuum provides the information on the 
light-meson spectrum, in particular the Goldstone boson-pion. However, it is difficult 
to solve directly the original theory based on the strongly correlated chiral-noninvariant 
vacuum. The general strategy is to perform a Bogoliubov-Valatin transformation and 
obtain a simple chiral-invariant vacuum at the expence of generating more complicated 
chiral-noninvariant quark interaction. Introducing concept of a quasiparticle, the quark 
mass becomes an unknown parameter, which can be determined by solving the gap 
equation. This program has been completed by Le Yaouanc et al. 0 and Adler and 
Davis [Q for the Goulomb gauge QGD Hamiltonian, without invoking explicit gluon 
degrees of freedom. As a result, the nonzero dynamical quark mass, m{k) ^ 0 , and the 
chiral condensates, 7 ^ 0 , as well as the tt — p mass splitting have been obtained. 

However, one encountered a ultraviolet divergent problem in both gap and Bethe-Salpeter 
equations. Depending on the formulation, one may also face the infrared problem. In this 
work, we include the dynamic gluon degrees of freedom utilizing the flow equations and 
obtain the renormalized gap and Bethe-Salpeter equations that are finite both in UV and 
IR regions. Introducing quasiparticles as elementary degrees of freedom, one may provide 
a possible connection between the BGS (pairing) model and the GQM. In particular, we 
investigate what is responsible for the n — p mass splitting. The spontaneous chiral 
symmetry breaking is nonperturbative in nature and manifests for small (zero) quark 
masses, while the hyperfine interaction is the leading order of the perturbation theory 
and describes nonrelativistic dynamics of large quark masses as dictated by the GQM. 

Our starting point in this approach is the QGD which has a complex chiral noninvari¬ 
ant vacuum and a strong chiral-invariant interaction between light particles. It is however 
an extremely difficult problem to solve the QGD at hand, because it is a strongly coupled 
theory with a nonfixed number of particles. Instead, we notice that the GQM provides 
reliable predictions in hadron physics and the relatively heavy constituents are coupled 
weakly by the chiral-noninvariant interactions. In the GQM, the vacuum is simple and 
invariant under the chiral transformations. It is an open question whether the weak cou¬ 
pling GQM is dual to the strong coupling QGD. Perhaps, an ultimate goal in solving 
QGD is to find such a dual transformation that converts the strong interacting QGD 
to a weak interacting effective theory. The many body approach suggests the BGS type 
Bogoliubov-Valatin (BV) transformation to fill the gap between QGD and GQM. We 
therefore adopt that the BV transformed QGD that has a chiral-invariant vacuum, but 
contains chiral-noninvariant interactions. The interactions can be decomposed into the 
strong BGS interactions, that reflect the chiral symmetry breaking and the residual weak 
interactions. Namely, the QGD Hamiltonian may be represented in the BV transformed 
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massive basis as 


HqcD — Hq + — {Ho + Hphen) + {Hqqd — Hphen) = Hpfp + Hpp , (1) 

where Hq is a free Hamiltonian, Hq^j^ contains all QCD interactions, and Hp^en is 
a phenomenological QCD inspired interaction Hamiltonian. We assume that after BV 
transformation, i.e. in the massive quasiparticle basis, one can choose Hphen in a way 
that it absorbes all stong interacting part of Therefore, {Hq + Hphen) = Hpfp is 

a strong interacting part which corresponds to the zero’th order BCS approximation, 
and {Hqqp — Hphen) = Hpp is a. weak residual interaction part which can be treated 
perturbatively. 

In this work, we use the Coulomb gauge QCD Hamiltonian, as a starting point. As 
discussed above, we then reorganize the Hamiltonian in such a way that the nonpertur- 
bative part includes the free Hamiltonian and the instantaneous interactions, summing 
the Coulomb and linear conhning potentials, Hpfp = Hq + Vinst, and the perturbative 
part corresponds to the dynamical interactions with propagating perturbative gluons, 
HpT = Vdyn- 

The success of the BCS (pairing) model 0 may be due to the explicit nonperturbative 
features such as dynamical chiral symmetry breaking and massive quasiparticle modes. 
Our aim here is to include the dynamical interactions perturbatively in the framework of 
BCS (pairing) model. Previously in Refs. P-0] the many-body relativistic Hamiltonian 
which contains the chiral-noninvariant interactions in a quasiparticle basis was solved 
directly for the lowest Fock sectors. In this way, only the particle-number-conserving in¬ 
teractions were taken into account. In principle, such a model 0 incorporates an exten¬ 
sive Fock space, because of the residual dynamical interactions which change the particle 
number to describe transition amplitudes between different Fock sectors. However, these 
interactions invariant under chiral symmetry were neglected since they were assumed to 
be weak. We now use the method of flow equations and include these supposedly weak 
interactions. 

The flow equations perform a sequence of unitary transformations to generate an 
effective Hamiltonian which is block-diagonal in the particle number space. This dehnes 
a resulting quasiparticle basis obtained both by the Bogoluibov-Valatin transformation 
and the flow equations, which includes the effects from the dynamical interactions as well 
as the dynamical chiral symmetry breaking of vacuum. Since the dynamical interactions 
mix low and high Fock states, they are responsible for the high energy region. Including 
these interactions is a vital step to reproduce a correct ultraviolet behavior, since it 
accounts for the dynamical propagation of gluons. This is missing when only the static 
chiral noninvariant interaction is taken into account. The presence of both static and 
dynamic interactions in an effective Hamiltonian makes it possible to compare our results 
with the covariant calculations. 

The paper is organized as followes. Section ^ concentrates on the formulation of 
an effective block-diagonal Hamiltonian in the Coulomb gauge QCD. In Section ^ this 
Hamiltonian is applied to one- and two- body sectors, and analytic formulations of the 
gap and Tamm-Dancoff/Random Phase Approximation equations are presented. In Sec¬ 
tion the gap and TDA/RPA equations are solved numerically and the results are 
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summarized. The summary and conclusions follow in Section |^. The details of the QCD 
motivated Hamiltonian, the second order flow equations, the gap equation with a double 
normal ordering, the analytic proof of nonzero pion mass in TDA/RPA and the RPA for 
the S and D wave p mesons are presented in Appendices 0,100 and ^ respectively. 


2 Effective QCD Hamiltonian in Coulomb gauge 

As mentioned in the introduction, we proceed with an effective Coulomb gauge QCD 
Hamiltonian given by 


H — Hq + Vinst + Vdyn — Hj^p + Hpp , ( 2 ) 

where Hq is the free Hamiltonian, Vinst is the instantaneous interaction describing static 
properties, and Vdyn is the dynamical interaction involving the gluon propagation. It 
may be considered in this Hamiltonian that the gluon field is decomposed into the non- 
perturbative, ^ and perturbative, components, i.e. = A^^ + Once 

the nonperturbative gluon configurations are averaged over after fixing the gauge and 
performing the BV transformation, the nonperturbative component may give rise to the 
strong instantaneous interactions both in quark P] and gluon sectors. The non¬ 
perturbative gluon field is then absorbed completely into Vinst and the nonperturbative 
gluon condensate (Hq) = (mixing terms with both A-^^ and A^'^ give either 

zero or higher powers in helds). On the other hand, the perturbative gluon field is in¬ 
cluded in the free Hamiltonian and the dynamical interactions. Both of them describe the 
propagating dynamical gluons. We combine the free Hamiltonian and the instantaneous 
interactions and build a nonperturbative part of the Hamiltonian, FIq -|- Vinst = H^p, 
which gives the zero’th order approximation. The dynamical interactions are included in 
the perturbative Hamiltonian, Vdyn = Hpp, and are treated perturbatively. 

In the Coulomb gauge, the free Hamiltonian is then given by 

Hq = J dxijj{x) {—ij ■'V + m)ip{x) 

+ Trjdx (U\x) + Bl{x)) , (3) 

where the non-abelian magnetic field is B = Bi = VjAk — VkAj + g[Aj,Ak], and its 
abelian part is represented by Ba- The degrees of freedom are the transverse perturbative 
gluon field A = A°‘T°‘ (A = A^^), its conjugate momentum 11, and the quark field in 
the Coulomb gauge. Motivated by the quark model phenomenology [], we introduce a 
conhning potential to account for the nonperturbative physics, and add it to the existing 
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Coulomb potential in the instantaneous interaction (only for the quark component) Q 

Vinst = -^Jdxdyij{x)'yoTy{x)VL+c{\x-y\)'ilj{y)'yoTy{y), ( 5 ) 

where the kernel is a sum of linear and Coulomb potentials dehned by 

CfVL+cir) = ar-Cf-^, ( 6 ) 

with the string tension a = O.lSGeld^ in accordance with lattice predictions, 

and the fundamental Casimir operator Cf = (A(? — l)/2iVc = 4/3 for A'/ = 3. The fourier 

transform of the instantaneous potential, Eq. (j^), gives 

n n h C/dTra, Sttu 

CfVL+c{\q\) = - ( 7 ) 

The dynamical interaction includes the minimal quark-gluon coupling, Vgg, and the non- 
abelian three- and four-gluon interactions, Vgg, i.e. Vdyn = Vgg + Vgg, where 

Vgg = -g J dxi>{x)-f ■ A{x)'iIj{x) 

Vgg = Trjdx (b\x) - B\{x)) , (8) 

with the perturbative gluon held A = and the corresponding perturbative compo¬ 
nent of the magnetic held. In what follows, we focus on the hrst term; the quark-gluon 
coupling. The goal is to use how equations to scale the dynamical interactions to lower 
energies by eliminating the quark-gluon coupling and to generate an ehective quark inter¬ 
action which can be diagonalized nonperturbatively for bound states when the generated 
terms are added to the instantaneous phenomenological interaction. In this process the 
renormalization will also be achieved by introducing second order 0{g‘^) canonical coun¬ 
terterms which remove the UV divergences from equations for physical quantities. Since 
the dynamical interactions represent the canonical Coulomb gauge QCD interactions 
proportional to the coupling constant, we assume that they are weak and the entire 
procedure may be conducted perturbatively. There is no over representation or double 
counting, because the phenomenological interaction determine the infrared (IR) behavior 
while the dynamical interactions free from IR divergences describe the UV part. As shown 
below, the dynamical part gives rise to the hyperhne interaction cx.T°‘'ip){gp'^cxT^'ip) in 
the quark sector, which should be added to the instantaneous terms in order to avoid 
the UV divergences and noncanonical renormalization in the gap and Bethe-Salpeter 
equations. 

^Generally, the nonabelian Coulomb interaction has both quark and gluon components 

Viu8t = -\j dxdyp°‘{x)V°-^{x,y)p'{y), ( 4 ) 

since the density p in the non-abelian Coulomb interaction is the full QCD color charge, given as 
p°‘{x) = 'ip(x)joT°‘'4’ix) -|- f°‘^‘^A^{x) ■ n'^(a:). 
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Along with the second quantization we choose basis states and expand the quark and 
gluon field operators in normal modes. Because of the phenomenological interaction in 
the quark sector given by Eq. (H), the trivial perturbative vacuum |0) may not be the 
minimum ground state for our Hamiltonian given by Eq. (§). Thus, we introduce a trial 
nonperturbative vacuum state |r2) containing quark condensate which can be determined 
variationally. This is analogous to the BCS type vacuum studies 0, 0. We refer the 
state |r2) as the BCS vacuum. The Fock space is constructed from this vacuum using 
quasiparticle operators and S which appear in the field expansions given by 

= / j^^[u{k,s)h{k,s)+v{-k,s)d){-k,s)]e^^^ 

aJ ^2uj{k) 

n(x) = -* E/ “) - “)!<=*". (9) 

where b\Q) = d\Q) = 0. Note that, since the phenomenological interaction appears only in 
quark sector, the ’gluon part’ of vacuum state is considered as trivial, i.e. a|r2) = a|0) = 0. 
Hence, the gluon Fock space represents eigenstates of free gluon Hamiltonian with single 
particle energy (jj{k) = |fc|. In Eq. (^, the quark operators are given in the helicity 
basis and all descrete numbers (helicity, color, and flavor for the quarks and color for 
the gluons) are collectively denoted as s and a, respectively. We use the spinors in the 
massive basis, i.e. the nonzero effective quark mass is explicitly included in the spinor; 

u(k,s) = A(*=) + ^ ) X. 

^ ^l + g (fc)x. \ 

\f2 V - s(fc)(cr • k)xs J 

v(-k, s) = ^E{k) + M{k) ^ 

1 ( s{k){a ■ k){-ia 2 Xs) \ 

^2 ^1 + s{k){-ia2Xs) J 


where the sine and cosine of the Bogoliubov angle ^(k) denoted as sin($(fc)) = s(k) 
and cos(^(k)) = c(k), respectively, are given by 


s(k) 

E{k) 


, c{k) = 

^k‘^ + M2(fc) ’ ^M2(fc) 

+ M 2 (fc). 


( 11 ) 


Here, E{k) is a single-quark energy which we refer as a gap energy. The effective quark 
mass M{k) is kept as an unknown variational parameter in the calculations and found 
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from the gap equation by minimizing the ground state (vacuum) energy. This approach is 
equivalent to set up the Coulomb gauge pairing model [0, that involves a Bogoliubov- 
Valatin transformation to a vacuum containing a qq condensate and uses optimization 
principle to give an equation for the condensate wave fnnction |r2(fc)). Here Xs and rjs are 
the standard two-component Pauli spinors of a particle and an antiparticle, respectively, 
and for an antiparticle t]_s = -io' 2 Xs, he. xi /2 = (1, 0), X-i /2 = (0,1) and t]_i /2 = (0,1), 
Vi /2 = 0)- With the definition given by Eq. ([T0|), the spinors satisfy the nonrelativis- 

tic normalization and orthogonality relations; {k, s)u{k, s) = v\—k,s)v{—k,s) = 1 
and u'^{k,s)v{—k,s) = v'^{—k, s)u{k, s) = 0. Canonical (anti)commutation relations are 

{b{k, s), b^{k', s')} = {d{—k, s), d^{—k', s')} = {27r)^S{k — k')6s^s' 

[ai{k, a), a](fc', a')] = (27r)^d(k - k')Dij(k)da,a', (12) 

where the gluon operators a = ai{k)°‘ = X]a=i ,2 ^(fc, A)a“(fc, A) are transverse, i.e. k ■ 
a“(fc) = k ■ a“l(fc) = 0, and Dij{k) is a polarization snm 

Ai(fc) = c(fc, X)€j{k, A) = 6ij - kikj , (13) 

A=l,2 

with nnit vector component ki = ki/\k\ and ki ■ Dij{k) = 0. 

The complete QCD motivated Hamiltonian, Eq. (|^), can now be expressed in second 
qnantized form, using Eq. (j^), and normal ordered with respect to the trial vacnnm |r2). In 
addition to the two-body interactions, normal ordering leads to the condensate (vacunm 
expectations) and one-body operators which should be regulated in UV region with the 
cut-off A S> ^QCD since they include the helds at one point. The hnal expressions are 
rather complicate as summarized in Appendix ^ 

Since onr goal is to obtain effective dynamical interactions in the qnark sector |, we set 
np the flow eqnation scheme for the Hamiltonian Eq. (^, ignoring pnre glnon non-abelian 
contribntion Vgg, i.e. for H = Hq + Vinst + Vqg. We restrict the nonpertnrbative part of 
the Hamiltonian to the diagonal (particle nnmber conserving) sector, since H^p is the 
zeroth approximation that should not be eliminated by flow equations. The pertnrbative 
Hamiltonian which includes the quark-gluon conpling and mixes different Fock sectors 
is attributed to the rest (particle number changing) sector; 

Hd = = Hq + Vinst 

Hr = V„. ( 14 ) 

Since the quark-glnon conpling is hrst order, 0(5'), we eliminate it using flow equations 
pertnrbatively. All terms are detailed in Appendix ^ along with the free Hamiltonian 
operator containing the kinetic K and the condensate O terms in the zero’th and second 

^ In this work only perturbative mixing between different Fock sectors via Vqg is considered. In order 
to consider the nonpertnrbative mixing between quark and gluon sectors, one should include the gluon 
component of the instantaneous phenomenological interaction and the dynamical triple- and four-gluon 
terms. Nonpertnrbative glueball calculations in a pure gluon sector can be found in the first paper of 

El- 
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orders in coupling, + 0^°^ (the same for h'q^). The details of Wegner’s flow 

equations in QCD were presented in the previous papers . 


Having identihed the diagonal and rest parts of the Hamiltonian, we now construct 
the generator, r], and solve flow equations for the hrst two leading orders. The leading 
order 0{g) flow equation reads 


dl 






( 15 ) 


where the kinetic and quark-gluon vertex terms are given by (see also Eqs. dm and 
(|164|) in Appendix 0) ; 


= E/ + 4(k)d,{k)] + -^i^(k)af{k)a‘{k ), 


dk 


(27r 


(27r)3 


(16) 


and 


''"-.S,/(. 4 iSi.) 

go{ki, k2, k3;l)ds^{-ki)T%s2{k2)-j^^^^vl {-ki)aiUs2{k2){27if6^^\ki - k2 - k^) 

^J2uj{ks) 

+ gi{ki,k2,k3;l)[bl {ki)T%s2ik2)-i^^^ul {ki)aiUs^{k2){27i)^6^^\ki - fcs - k^) 

- d\^{-k,)T<^ds,{-k2)^^M=vl{-k2)a,VsA-ki){2n)H^^\k2 -k,- k,)] 


2 a; (^3 


gr{ki, fca, k^; 1) 4,(-fc2)T°6^3(fc3)n]^(-fc2)ain^3(fc3)(27r)^^^^^(fc3 - fci - fca 

'2a;(fci) 


H.c. 


( 17 ) 


respectively. Here, E{k) = ^Jk"^ + M‘^{k) (Eq. ([TI|) ) anda;(A;) = k. By implementing flow 
equations, the effective coupling constants are generated (see below) as go{ki, k 2 , k^; 1), 
gi{ki, ^ 2 , k^] 1) and gii{ki, ^ 2 , k^] 1), which are the functions of all three momenta corre¬ 
sponding to a given Fock sector (Note here that the different Fock sector operators flow 
differently in energy-momentum as indicated by 0, 1 and 1' subscripts here and below) 
and depend upon the flow parameter /. From Eq. (|T7|) the leading order generator is 


= — 


E /(n^ 


si,S2,a \n 


?7o(fci,fc2,fc3;/)4^(-fci)T°5^,(fc2) {-ki)aiUs2{k2){27iY5^^\ki ~ k2 - k^) 

2u{k3) 































- <(-fcOTX.(-fc2) f^^^\ vl{-k2)a,VsA-ki){2'K)H^^\k^ -k,- k,)] 


s/2u{k3 


r7i,(fci,fc2,fc3;/)4=^^^2(-fc2)TX(fc3)^^l,(-fc2)«*M.3(fc3)(27r)3<5(3)(fc3 _ k, - k^) 

l2u{ki) 


- H.c. 


( 18 ) 


where 


r]i{ki,k2,k3-,l) = Di{ki,k2,k3)gi{ki,k2,k3-,l), ( 19 ) 

for i = 0,1,1', with energy terms (i.e. energy denominators in the old fashioned pertur¬ 
bation theory) 


Do{ki,k2,k3) — —(E(ki) + E(k2)+uj{k3)) 

Di{ki,k2,k3) = E{ki) - E{k2) - u{k3) 

Dy{ki,k2,k3) = u;{ki) - E{k2) - E{k3). ( 20 ) 

The solution of the flow equations, Eq. (^), for the effective coupling constants is given 
by 


9S) = f/( 0 )exp(-T) 2 /) ^ ( 21 ) 

which eliminates the quark-gluon coupling, Eq. ( 0 ). for Z —>• cx), as anticipated. Corre¬ 
spondingly, the new operators in the particle number conserving sectors are generated 
from the second order 0{g‘^) flow equation 

( 22 ) 

which contribute to the effective block-diagonal Hamiltonian. One can always eliminate 
the particle number changing off-diagonal terms, H^‘^\ appearing in second order by 
choosing the generator in addition to Solving the flow equations, 

Eq. (|l^ and Eq. (^2]), the block-diagonal effective Hamiltonian, hfe//, renormalized to 
the second order is obtained (See details in Appendix P) as 

He//(A) = Ho{k) + V,nst{^) + VgU^)+6XcT{^), (23) 

where Vgen includes new operators generated via perturbative elimination of the dynam¬ 
ical interaction, the quark-gluon coupling Vqg, present in the original Hamiltonian, Eq. 
(H). The cut-off dependence on A is introduced by regulating operator products which ap¬ 
pear at one point and contribute UV divergences in the momentum space loop integrals. 
Regularization results in a large sensitivity of the effective Hamiltonian to parameter A 
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which is set at the UV scale. In order to eliminate this A dependence, a renormaliza¬ 
tion procedure is used. The renormalization is achieved through second order by adding 
the counterterm, 6Xct, such that in the limit A —>■ cx) the effective Hamiltonian, Eq. 

and hence equations derived using H^ff are UV finite and do not depend on the 
cut-off parameter A. Therefore we can omit the cut-off notation in H^ff, Eq. (^3]); i.e. 

= H,„. As shown below, the cut-off sensitivity is very weak in the perturbative 
renormalization. 

In what follows we focus on the quark sector and consider only the quark renormal¬ 
ization which goes through the second order and includes quark self-energy calculation. 
In the Hartree approximation the Dyson equations for the instantaneous and generated 
self-energies are (See details in Appendix P) 




where, using Eqs. (§) and (|T0|) , the equal time 3-d quark propagator, S^‘^\ is related to 
the Feynman propagator, 

^ ^ ~ ’’’' ~ + ■>'' kc(k))ag'j 


and is given by 


S®(/!) = J^iSf‘Hk,k„) = ^{s(k)-Tkc(k)) . 
In Eq. the potential functions are given by 


CfVL+cik, q) 
CfW{k,q) 


1 Cfg'^ Aira 

2 {k — qy (fc — qY 

__ 

u;{k - q){E{q)+u;{k - q)) ' 


(26) 


(27) 


Note that this form of generated potential is valid for a large loop momentum q and 
sufficient for the renormalization. In general W is more complicate (For details, see 
Appendices P and P). The complete proper self-energy reads 


S 


-‘inst 


+ E 


gen ■ 


( 28 ) 
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Making a non-relativistic ansatz for the self-energy, E can be represented as 


S = S(fc) =mkl(fc)+7-fc5(fc), (29) 

where A and B are some scalar functions. In order to hnd A we apply trace, Tr, to the 
l.h.s. and r.h.s. of Eq. (p^). Multiplying the hh.s. and r.h.s. of Eq. (p^ by 7 ■ fc and 
applying trace after that, one can hnd B. The following self-energies are obtained 

^inst(k) = J j^^CfVL+c{k,q)s{q)e-^^/^^ 

+ T ■ y j0^CfVL+c{k, q)c{q)k ■ 

S,en(fc) = j ^^CfW{k,q)s{q)e-^^^/^^ 

+ T ■ q)c{q)k ■ iq ■ , (30) 

where I = k — q, and exponents are the regulating functions (Appendix P). We calculate 
the divergences associated with the self-energy S (Eq. (pQ])) given by 

= mA^^^ + j ■ kB^^\ (31) 

and dehne the quark mass and wave function corrections as 

6m = mA'^'^ 

Z -I = B '^^^, (32) 

respectively. | Separating the divergences, one has 

S£,(fc) = rnj ^l\mC,VL+c(k,q)e-'’"/''" 

+ ' '“I / 7^ 11“ «)* ■ 

since for the large momenta, q —> 00 , s{q) —> m/q and c{q) 1 with the current 
quark mass m. The integrals following the structures Im and 7 ■ k are identihed as 

^ In general, 

5m = 

Z-1 = (dE^“/(7^fc'')) , (33) 

where fh is the renormalization point (the renormalized mass). This definition is consistent with the one 
given above. 
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j^div respectively. In Eq. (^), includes the instantaneous and generated 

divergent contributions, while only a Coulomb potential in the in¬ 

stantaneous interaction has UV divergent behavior. As shown below, the divergent terms 
j^div j^div hence both corrections 5m and Z reduce to constants. This means 
that the effective Hamiltonian Eq. (P3|) can be renormalized canonically by introducing 
a momentum-independent counterterm, SXct, which corresponds to a local operator. 
The divergent self-energy operator, Eq. ( 0 ). appears as a correction to the quark free 
Hamiltonian, Eq. 

= J dxip{x)[{m + 5m) — Z■'V]^jJ{x), (35) 


that justihes the above dehnitions of the mass and wave function corrections, Eq. (|3^) . 
Calculating integrals in Eq. (^), we hnd 


^div 

■^inst 


^div 

^gen 



4 In A j -|- 7 ■ fc 
2 In A ) -|- 7 ■ fc 


iFi"*) 


(36) 


where, as expected, the leading A^ divergence does not appear in the quark sector. This 
results in 


5m = In A 

(dvr) 2 

Z = 1, 


(37) 


for the sum of instantaneous and generated terms. Note that the combined instantaneous 
and dynamical terms together do not require the wave function renormalization, however 
each term alone requires this renormalization. Therefore the only necessary counterterm 
is the quark mass counterterm 

5Xct{A) = Mct{A) j dx^{x)'ip{x ), (38) 


which absorbs the UV divergences in the quark sector of i7e// as A —> cx). Here, 


Mct{J5) = —5m = 


Cfl 

(dvr)" 


6m In A. 


(39) 


The mass counterterm is proportional to the bare quark mass, m, and thus vanishes in 
the chiral limit m —> 0. This means that, provided both instantaneous and dynamical 
terms are included, the quark sector of the effective Hamiltonian, Eq. (^31) , is UV hnite 
in the chiral limit and does not require any renormalization at all. As shown above, Eq. 


(|36[) , this happens due to a complete cancelation of divergent wave function corrections 
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from the instantaneous and generated terms. The renormalized free quark Hamiltonian 
reads 


= Hq + SXct{-^) = J dxi/j{x)m{A)i/j{x) 


m(A) = m + Mct{^) = m ( 1 — 


Cfl 

(dvr)" 


6 In A 


(40) 


which insures that the effective Hamiltonian is UV hnite in the quark sector. 

Our results are in accordance with the canonical renormalization of the Coulomb 
gauge QCD Hamiltonian [|T^], Hqcd, since both Hamiltonians, Heff and Hqcd, have 


the same UV behavior. Extensive discussion in the literature is devoted to a correct 
renormalization procedure of an effective Hamiltonian based on the Coulomb gauge, and 
the corresponding formulation of UV hnite equations 0. The authors in Ref. 0 claim 
that in the chiral limit one is unable to avoid UV divergences and therefore suggested to 
introduce a non-canonical momentum-dependent counterterm to absorb them. In related 
studies the authors adopt a special prescription of double normal ordering with respect 
to the perturbative, |0), and the non-perturbative, jH), vacua in order to cancel the 
divergent term (see Appendix |^. However, as shown in Appendix P, this prescription 
fails for some potentials in the IR region. In all of these cases, the divergence is caused 
by the term P] 


which appears from the instantaneous interaction in the chiral limit (see Eq. (|3^). As 
shown above, the dynamical interaction which is missing in the mentioned works i, 0 
has to be included in the Hamiltonian in order to cancel the instantaneous divergent 
contribution. It may be analogous to the standard time-ordered perturbation theory, 
where all time-ordered diagrams in a given order should be added to obtain a correct 
CO variant result. 

We summarize the matrix elements of Hg/f, Eq. (^), in the quark sectors of interest 
(up to two quark states) 

= 0--(A) + 0,„,i(A)+0,en(A) 

= K^^^{A) + J:^nst{^) + ^geniA) 

(2|-f^e//|2) = Vinst + Vgem (42) 

where IH) is a shorthand notation for the zero-quark sector (also vacuum state), |1) is the 
single-quark sector, etc., and all other terms are specified below. One should distinguish 
between two types of i^e// terms. The hrst one arises from the normal ordering of the 
original Hamiltonian, Eq. (^: the instantaneous interaction with linear plus Coulomb 
potentials, labeled by inst. This leads to the self-energy operator T^inst in the one-body 
sector, and the condensate term Oinst in the zero-body sector. The energy of the quark 
ground state O = Oq comes from the normal ordering of the free quark Hamiltonian 
Hq with respect to the vacuum |r2). The second type of terms are dynamical operators 
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generated by flow equations and labeled by gen. In Eq. (^21) the renormalized condensate 
O'’®"' and the kinetic terms are given by (Appendix 


0’'®"(A) = O + 5X°“"(A) = -4iV,y I ^[kc{k) + m(A)s(fc) ] 

= + rn{A)s{k))[bl{k)b,{k) + dl{k)ds{k)] 

+ (fe(fe)-m(A)c(fe))|6l(fe)(it(-k)+(i,(-k)6,(fe)]] , 


(43) 


where m(A) is deflned in Eq. (^). Here, K = Kq and O = Oq are deflned by Eqs. ([149| ) 
and (|154| ), respectively, in the Appendix ^ and is the mass counterterm given by 

Eq. (^) in the zero-body sector (analogous for 5X^^'^^). The mass counterterms 5Xct 
cancel the leading In A behavior of the radiative corrections to the vacuum and kinetic 
terms. In Eq. (|4^) the corrections to O'”®" and regulated by the exponential cut-off 

function, include the condensate terms (Appendix 


Oinst{A) 


Ogen{A) 




1 - s{k)s{q) 


c{k)c{q)k ■ q 


e 


1 -|- s{k)s{q) -|- c{k)c{q)k Iq 1 


(g+fc)VA^ ^44) 

^-(g+fc+OVA" 


and the polarization operators (Appendix P) 


Sinst(A) — 


dkdq 

(^ 


CfVL+cik,q) s{k)s{q) + c{k)c{q)k ■ q 


X [bl{k)bs{k) + dl{-k)ds{-k)] 

r dkdq 

^ sJ (27r)6^^ 

X [6j(fc)(ij(-k) -F dsi-'k)bs{k)] 


CfVL+c{k, q) -c{k)s{q) + s{k)c{q)k ■ q 


E,,„(A) = y: 


dkdq 

(^ 


CfW{k, q) s{k)s{q) -|- c{k)c{q)k Iq l 




,-qVX 


.-V/A 2 


X [bl{k)bs{k) + d\{-k)ds{-k)\ 

^:CfW{k,q) —c{k)s{q) + s{k)c{q)k ■ iq ■ i 


E 


(27r) 




X [blik)dl{-k) + dli-kMk)], 


(45) 


where the potential functions lA+c and W are deflned in Eq. (0. The effective quark 
interaction includes the two interactions, Vinst + Vgen, that define the effective Hamilto¬ 
nian Eq. (1^) in the two-body sector, (2|ife//|2). We consider only two-quark interactions 
which contribute to a meson bound state equation in Tamm-Dancoff (TDA) and Ran¬ 
dom Phase (RPA) approximations. In the c.m. frame the instantaneous and generated 
interactions contributing to TDA (X component of the RPA wave function) are (see 


14 



























Appendix |^) 




X 

+ 


^gen 


X 

+ 


(4, 

{ul{q)ua{k)){vl{-k)v^{-q)):bl{q)T%a{k)df3{-k)T‘^dl{-q): 

{ul{k)us{q)){vl{-q)vp{-k)):bl{k)T%s{q)d^{-q)T‘'dl{-k): 

S /S2»'.(fc.<;)Aj(fe-9) 

{ul{q)aiUa{k)){vl{-k)ajV^{-q)):bl{q)T%aik)dfs{-k)T'"d\{-q): 
{ul{k)aiUs{q)){v\{-q)ajVp{-k)): bl{k)T%siq)d^{-q)T'^dl{-k): . 


In RPA {Y component of the RPA wave function), they are given by 


K; 


inst 


^gen 


■ SJ 

X {vl{-q)Uo,{k)){vl{-k)us{q))-. d^{-q)T%^{k)d0{-k)T%s{q)- 

+ {ul{q)v/3{-k)){ul{k)v^{-q)): bl{q)T'^dl{-k)bi{k)T<^dl{-q): 
dkdq^ 


= E 

a.j35^ 


(2vr) 


, 2 W 2 {k,q)Dij{k - q) 


X 


{v\{-q)aiUc,{k)){vl{-k)ajU5{q))\d^{-q)T%a{k)di3{-k)T%s{q)- 

where potential functions are given 


(47) 


CfWr{k,q) 

CfW2{k,q) 


1 __ 

2uj‘^{k-q) + {E{k)-E{q)y 

1 _ Cfl _ 

2u;^{k-q) + {E{k) + E{q))^ ' 


(48) 


and Vl+c is dehned in Eq. (P7|). In the next section we utilize the obtained effective 
Hamiltonian, with matrix elements given by Eq. (^^, to derive and solve the quark gap 
and meson bound state equations. 


3 Sector solution of the effective Hamiltonian: gap 
and bound—state equations 

Now that we have eliminated the quark-gluon coupling, which mixes different quark sec¬ 
tors, and obtained the effective quark Hamiltonian, H^ff, valid up to the second order 
in the coupling constant, we can nonperturbatively diagonalize each sector Hamiltonian, 
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Eq. (|^). Further, because -ffe// is also renormalized, the equations for physical observ¬ 
ables are free from the UV divergences. In subsection p.l| we hrst investigate the quark 
vacuum by formulating the quark gap equation and also calculate the quark condensate. 


Then we address the meson spectrum in subsection |3.2| . Numerical solutions of the gap 
equation and the bound state equations are discussed in the next section (Section HI). 


3.1 Gap equation 

The gap equation allows the determination of a nontrivial vacuum with quark conden¬ 
sates and propagating quasiparticles (or quarks with a dynamical mass). There are several 
ways to obtain this equation, the most common based upon a variational principle to 
minimize the vacuum (ground state) energy. The variational parameter is the angle of 
transformation from undressed to dressed particle (quasiparticle) operators, <h(fc), which 
dehnes a quasiparticle basis, Eq. (0, with a dynamical quark mass M{k). Therefore, 
minimizing the vacuum energy of the effective Hamiltonian, i.e. 


Siik) 


= 0 , 


(49) 


generates the gap equation for the unknown <h(fc) or M{k). Using Eq. (^) for the 
condensate terms 


(Or 


6^{k) 

the following gap equation is obtained: 

dq 


Oi 


inst 


o 


gen) 


= 0 , 


(50) 


ks{k) — m{A)c{k) = 

f dq 


(27r)= 


CfVL+c{k,q) c{k)s{q) - s{k)c{q)k ■ q 




(27r)= 


CfW (fc, q) c{k)s{q) — s{k)c{q)k ■ Iq 1 


.-V/A 2 


(51) 

, and m(A) 


where I = k — q, potential functions Vl+c and W are given in Eq. 
includes the mass counterterm (Eq. (^)) dehned by Eq. (^^. 

The gap equation can also be obtained by demanding that the effective Hamiltonian 
should not contain off-diagonal one-body terms of the type bd and b^dK This means that 
the BCS vacuum |f2) is stable against quasiparticle pair creation. Therefore, the operator 


E 


dk 


s ^ (27r)3 

of Heff vanishes by choosing 


F(fe,®)|6;(fc)4(-fe) +<i,(-fc)6.(fe)] 


F(k,<S>) = 0. 


(52) 


(53) 


Imposing this condition on the nondiagonal matrix elements, labeled nd, of the effective 
Hamiltonian in the single quark sector (l|ife//|l) given by Eqs. (|^ and (^), i.e. 


+ j:]^/j\nondiag. 


T/'Tid I ^nd 

^ren + ^4 


I ynd _ Q 
inst ' ^gen ^ ’ 


(54) 
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yields the same gap equation as above. 

The alternative way to obtain the quark gap equation without specifying matrix 
elements of hfe// is to use the Dyson equation for the self-energy operator S, Eq. 


With the propagator given by Eq. the instantaneous and generated self-energies 
have been found in section 0, Eq. (^). Therefore A and B functions, dehned in Eq. 


are 


kB(k) = j -i0j^CjVL+c(k,q)c(q)k-qe~'’^''^' + j -^^CfW{k,q)c(q)k-tq . 


On the other hand, using the general expression for the self-energy S, Eq. (PU]), the 4-d 
dressed Feynman propagator is 


S^^\k,ko) 


D(fc) 


1 

7 o^o — 'y ■ k — m — E(fc) 

'yoko - ^ ■ k{l + B{k)) + m{l + A{k)) _ R+{k) R_{k) 

ki-n‘^{k) “ ko - n{k) ^ ko + n{k) 

\jk'^{l -l- B{k)Y + m^(l -|- A{k)y , (56) 


with the residues R±{k) given by 


R±{k) 



m(l -|- A{k)) 
D(fc) 


T -7 ■ k 
2 ’ 


1 + B{k) 

n{k) 


(57) 


The integration over ko, Eq. ([25|), with the residues R±{k) gives the 3-d equal-time 
propagator 


S^^\k) 


f dko ^ _ m(l + A{k)) - 7 ■ fc(l + B{k)) 

J “ Wk) 


and comparing Eq. (|5^) with Eq. (p^ for we see that 


s{k) 


m{l + A{k)) 

D(fc) 


, c{k) 


k{l + B{k)) 

n{k) 


(58) 


(59) 


where Q{k) is given in Eq. (|56|) . Essentially Eq. (|5^ is the gap equation. Eliminating D 
in Eq. (^ 


s{k)k{l + B{k)) — c{k)m{l + A{k)) = 0, (60) 

and substituting A and B, Eq. (^), we obtain again the same quark gap equation, Eq. 

(ID- 

Due to the gap equation, Eq. (|MD , the single-quark operator is diagonal 

H + dl{k)ds{k )], (61) 
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and therefore it can be associated with an effective qnark energy, e{k), 


H, [ Hr + EI ^ Ki„ + Yin., + Kn = ■ (62) 

The vanishing non-diagonal part of this operator, Eq. (M), is given by Eq. (|60D , then its 
diagonal part is 


e{k) = c{k)k{l + B{k)) + s{k)m{l + A{k)), 
and substitnting A and B fnnctions, Eq. (pSD , gives 

dq 


e{k) = kc{k) + m{A)s{k) + 


(27,) 


;C’fVi,+cik, q) s(k)s(q) + c(k)c(q)k ■ q 


dq 

W)- 


■CfW{k, q) s{k)s{q) -|- c{k)c{q)k Iq l 




(63) 




(64) 


This equation can also be obtained using (l|hfe//|l), Eqs. (^2]) and (|^. It is convenient 
to represent e{k), using the gap equation, Eq. (|6^), k{l + B) = [c{k)/s{k)]m{l + A), as 


e{k) = 


m(l -|- A{k)) 


s{k) 


m 


dq 


s{k) J (27r) 


-CfVL+c{k,q) 


s{q) 

s{k) 




(27r 


s{q) 

s{k) 


(65) 


The effective energy e{k), Eq. (p3|), equals to the frequency fl{k), appearing as a pole in 
the propagator, Eq. (^), 


e{k) = fl{k). (66) 

In order to introduce e{k) the Bogoluibov-Valatin angle should satisfy the gap equation, 
Eq. (|59|) or Eq. (|60|) . Substituting Eq. (|^) into e{k), Eq. (|6^ , gives Eq. ([661) . However, 
one should distinguish between the effective energy e{k) and the gap energy E{k), Eq. 
O- Comparing Eq. m with Eq. (|5^) for the sine and cosine we hnd 

E{k) = eik)/{l + B{k)). (67) 


A natural question arises, what value can be considered physical, the effective quark 
energy e{k) = fl{k) or the mass gap M{k) (or related gap energy E{k)). Consider Erst 
the propagator for a free massive Dirac particle with mass M{k) 


1 _ joko - ^ ■ k + M{k) 

7o^o — 7 ■ — M{k) fcg — E‘^{k) 


R{k) 

ko-E{k) 


-|- [analytic at ko = E{H^S) 


with E{k) = + M‘^{k) and the residue R{k) at the positive frequency pole given by 


R{k) = 


7o + 


M{k) 

'W) 
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where we have used Eq. (|TI| ) for the sine and cosine. Let us compare this with the 
propagator in our model, Eqs. (|5^) and (^7|) 

E(fe) = 1^^) + atk„^a{k)], (70) 

where, using Eq. for the sine and cosine, the residue R+{k) at the positive frequency 
pole is given by 


R+{k) = - 


7o + 


m(l + A{k)) 


- \l ■ ^ [70 + s{k)] - ^7 ■ kc{k ), (71) 


which is the same as Eq. ([691). Using Eq. (^) and Eq. 


propagator Eq. (|70D can be written as 


70^0 - (7 ■ fc + M{k))Vt{k)/E{k) 


for the sine and cosine, the 


(72) 


Hence, though the two covariant 4-d Feynman propagators, Eqs. (^) and (|70D (or (0)), 
have different behavior with different poles, E{k) and Vt{k) = e{k), respectively, they 
have the same 3-d image S^^\k), given by Eq. (^ 6 [) , which appears as an equal-time 
quark propagator in the Coulomb gauge Hamiltonian, Eq. (H). Therefore, E{k) can be 
considered as a physical pole, associated with a quasiparticle having an effective mass 
M{k). This also justihes our dehnition of sine and cosine through M{k) and E{k), Eq. 
O- We show below that the frequency pole Q{k) is not well dehned, and thus can not 
represent a particle. 

We investigate the UV and IR behavior of e{k) = fl{k), Eq. (|64D, and M{k), given 
by the gap equation Eq. (|^). In the UV region both equations are hnite, since based 
on the renormalized Heff{A), the mass counterterm m(A) cancels exactly the only UV- 
divergence from A{k), Eq. (^), while B{k) does not contribute to the UV-divergence 
(as discussed in Section In the chiral limit e{k) and M{k) are dehned by the UV 
hnite equations even without renormalization (See Section ^). 

Infrared problems are caused by the linear potential in the instantaneous interaction, 
V^^cik, q), diverging as \k — q\~^ as k ^ q. Consider the limit k ^ q in the gap 


equation, Eq. (|^), 


dq 


J (27r)3 

Expanding up to the second order 


CfVL+cik,q) c{k)s{q) - s{k)c{q)k ■ q 


(73) 


s{q) = s{k) + k ■ ds'{k) + 0(<5^) 
c{q) = c{k) + k ■ Sc'{k) + 0(<5^) 
kq = 1 + 0(6“^), 


(74) 


where 6 = q — k and s'(fc) denotes the derivative of sine in k, we hnd for the term Eq. 

o 




k ■ S[c{k)s'{k) 


sikyik)] + ois^) 


(75) 
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which behaves, after angular averaging, as 


dq\k-q\ 0{{q-k)), 


(76) 


that converges. Thus, the mass gap M{k), Eq. (^), is well defined in the IR. In Appendix 
y we discuss the gap equation obtained by double normal ordering with respect to the 
perturbative |0) and nonperturbative |fl) vacuum states (See Eqs. (|207|) and (p08|) ) which 


was used in Ref. as well as the first paper in Ref. |^. One of the motivations for this 
prescription was to avoid UV divergences for the Coulomb potential in the chiral limit. 
However, the dangerous term reads 


dq 

(2vr)- 


■CfVL+c{k, q) c{k)s{q) - s{k){c{q) - l)k ■ q 


(77) 


which, using Eq. (0. behaves as fc —> q as 

f dq 


(2w) 


rC/U+cCfc, q)s(k) 


dq\k — q\ ^s(fc) , 


(78) 


that diverges. Thus, the double normal ordering gap equation is IR singular, and a mass 
gap M{k) does not exist for a confining potential. On the contrary, our gap equation, 
Eq. (^Ij), is well defined for the Coulomb plus linear potentials in the UV and IR, and 
provides a finite mass gap M{k) (or a finite gap energy E{k)). 

However, the effective quark energy e{k), Eq. 


dq 

W)- 


■CjVi,+c(k, q) s(k)s(q) + c(k)c(q)k ■ q 


(79) 


and the frequency pole ff(A!), Eq. (|6(^), behave for a confining potential as fc -e q as 

f dq 


(27r) 


rC/I4+c(fc, q) 


dq\k — q 


-4 


(80) 


that diverges. Hence, contrary to the findings in Ref. 0, the effective energy is not a 
physical observable and is clearly not IR finite. Instead the excitation energy e(k) — ^(0) 
is IR finite. Using Eq. (|64D and Vi+cik, q) —>■ Vij^cik — q), 

dq 


= I 7 KZ^^f^L+c{k,q) s(fc)s(q) + c(fc)c(q)fc • q +IR finite (81) 

-CfVL^ciq) s{k)s{k — q) + c{k)c{k — q)k ■ {k — q) +IR finite, 


{271 f 

dq 


J (2rr) 

the excitation energy is given by the IR finite formula 

dq 


e{k) — e(0) = 


J (27r 
IR finite 


-CfVL+c{q) s{k)s{k - q) - s(q) + c{k)c{k - q)k ■ {k — q) 


dq\q\ 0(q), 


(82) 


20 























as q —0. The same holds for fl{k) — f2(0). Using Eqs. (pQ]) and 

dq 


m(l + A(k)) 

n{k) = + 


s{k) 


(2vr; 


-CfVL+cik, q)s{q) / s{k) + infinite 


m + j -^^CfVL+c{q)s{k - q)]/s{k) + infinite, 


( 2 vr) 


( 83 ) 


the frequency pole difference is given by 

dq 

(27 

dq\q\~'^0{q^) , 


n(fc) - (1(0) = j ^C;V,+o(q) 


s{k - q) 
s{k) 


s[q) 


+ infinite 


(84) 


that converges as q —0. However, by shifting the pole, the Feynman propagator, Eq. 
(1^, still has the IR divergent structure 




R+(k) 


+ 


n{k) 


fco - H(0) - (H(A;) - H(0)) fco + H(0) + (H(fc) - H(0)) ’ 


(85) 


containing the IR divergent term f2(0) in the denominators but with the residues 
given by the IR hnite expressions, Eq. o, 

R±(k) = J l7o±s(fe)l T j7'Mfe) ■ (86) 

This structure is a reflection of the fact that as a result of conhnement, an inhnite amount 
of energy is required to create a single quasiparticle state from the vacuum. This means 
that a color singlet state does not exist and cannot be a physical state in our model. 
This is generally true in hadron physics. At the same time the energy gap, given by the 
gap equation, Eq. (|^ , does exist and corresponds to a physical quantity, dehning the 
gap between a vacuum |H) and the hadron scale. This interpretation makes possible to 
map our model on the constituent quark model by associating the energy gap with the 
elementary degrees of freedom, quasiparticles. Quasiparticles with effective mass M{k) 
correspond to the valence quarks. The dynamics of quasiparticles is described by the 
Feynman propagator, Eq. (|68[), which is IR hnite with a physical pole at M{k). In 3-d 
a pole of the equal-time propagator is given by the effective energy of a quasiparticle, 
E{k) = ^k‘^ + M^{k). 

Next we consider the quark condensate {Q\'ijj'if\Q) for a single quark havor. Using Eqs. 
(H) and (|l3), we obtain 

{nifjifin) = J = -2N^ J J^s{k). (87) 


This can also be evaluated in terms of the equal-time propagator S^^\k), 

_ 1 _ 1 _ p 

= A'=fe7{Si|-fe(o), v>o(o)| + -{</)j(o),^4o)}|sj) = - n^j 
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where the minus sign arises because l/2['0,'?/>] —>• Substituting Eq. (p6|) for S^^\k) 

into Eq. (|88|) , we get the quark condensate of Eq. (^7|). We regulate the quark condensate, 
Eq. (|87|), by subtracting the perturbative contribution 




-w 


-2W 


dk 

(27r)3 

r dk 


i27r] 


(Tr5®(fc) -Tr5®(fc)) 


m 




(89) 


where m is the bare quark mass. As |fc| —>■ cx) the mass gap M{k) —>■ m and the 
nonperturbative sine behaves as s{k) —>■ m/y/k"^ + Thus this subtraction improves 
the convergence of the quark condensate integral in the UV. 

For further investigations it is convenient to introduce the scalar, Us, and vector, cx^, 
parts of the equal-time quark propagator, in Eq. (|^), i.e. 

S'(3) ^ ^ ('9Q) 


and express all quantities of interest using them. In particular, in the chiral limit m 
one has from Eqs. (|^ and 


0 , 


<{k) 


Mo(fc) 

W)o 

Eo{k) 

klo{k) 


^ 1 

2n{k) 2 ^ ^ 

1 + B{k) _ lc(fc) 
2n{k) ~ 2~T~ 

mA{k) <Js{k) 
1 + B{k) (Ty{k) 


-ON [ ^Mk) 
"7 (27r)3 n{k) 

n{k) _ 1 

1 + B{k) 2cr^(fc) 



£(k) = ^m'^A'‘[k)+B[l + B(k)Y, 


(91) 


where A and B functions are dehned in Eq. (|5^) , and subscript 0 denotes the chiral 
limit case. When the scalar part mA ^ 0, the nonzero mass gap Mq = mAj (1 -|- B) and 
chiral condensate ('0V’)o ~ / dkmA/Q are generated, and there is a dynamical chiral 
symmetry breaking (m = 0 corresponds to no explicit chiral symmetry breaking). An 
obvious property following from Eq. (E3) 

trj(fc) + k'^aUk) = i (92) 


shows that low momentum behavior is governed by the scalar part of the propagator, ag- 
One can generalize the above expressions for a nonzero current quark mass m by 
substituting mA{k) m{l + A{k)) in Eq. (|9T|) , and regulating the condensate by 
subtracting the perturbative value. 
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3.2 Meson bound state equations: TDA and RPA 

As discussed in Section we map our effective model on the constituent quark model 
with the quasiparticles of mass M{k) playing the role of the valence quarks. In this Sec¬ 
tion, we represent mesons as bound-states consisting of quasiparticles and seek approxi¬ 
mate eigensolution of our effective Hamiltonian. We do not solve here the Bethe-Salpeter 
or Salpeter equations. Instead, we formulate bound state problem in the Tamm-Dancoff 
approximation (TDA) and subsequently in the Random Phase (RPA) approximation 
[p0| . We show below that the flavour octet family can be described well within these ap¬ 
proximations, separating the chiral symmetry and spin effects. However, flavour singlet 
states need more elaborate approach, which we discuss in a separate publication. 

In terms of the quasiparticle operators used in Section the TDA meson creation 
operator reads 

Ri = j . (93) 

which, acting on the vacuum, creates a meson with a wavefunction \'ipn) with the quantum 
number n, and annihilates into the nonperturbative vacuum, 

^nl^) = IV’n) 

= 0. (94) 

Eq. (1^) can be considered as a dehnition of the TDA vacuum, consisting of a condensate 
of only quasiparticle-quasihole pairs. The commutation relation of the meson operators 

= (95) 

leads to a normalization condition for the wave functions 

/ 7^ II {q)i^n (q) = N6nn' , (96) 

5-y 

where N is the normalization constant. Projecting the Schrodinger equation i7e//|V’n) = 
Enli^n) onto one-particle-one-hole truncated Fock sector, we get the TDA equation 

where the binding energy is dehned as — Eq with the vacuum energy Eq 

subtracted {Heff\fl) = Eo|D)) and in the r.h.s. of the TDA equation the operator 
picks up the wave function component 

{Ji|K,g(fc)4(-fc)]|Ji) = (98) 

The TDA is improved by extending the quark vacuum to contain also four quasiparticle 
condensates in addition to two quasiparticle condensates. Including up to four quasipar¬ 
ticle correlations beyond the BCS is known as the RPA approach. Generalization of the 
operator of Eq. (p3D containes meson creation and annihilation terms 

Qi = I . (99) 
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The RPA wavefunction and the RPA vacuum are 


Qnl^) = IV’n) 

Qn|f^) = 0 , ( 100 ) 

where, though the same notations were used as above, they should not be confused with 
the TDA wave function and TDA vacuum. From the meson commutation relation 




( 101 ) 


the following normalization condition for the wave function components, X and R, is 
obtained 


/ 


dq 


5'y 


= N6 


nn' 1 


( 102 ) 


with the normalization constant N. To derive the RPA equations of motion we calculate 
the commutators 


(Sl|lQ„,|ff.„,6j.(fe)4(-fc)]||Sl) = M^Xf(k) 
(Sl|ll3„,|/f,„,()„(fe)<is(-fc)]||Sl) = (103) 

which pick up the X and Y components in the r.h.s. of equations 

(Sl|[Q„,4(fc)4(-fc)]|li) = xfik) 

{a\\QnXi.k)d^{-k)]\a) = Yf(k). (104) 


The RPA system of equations, Eq. ( |103| ), reduces to the TDA, Eq. (p7|), by putting 
R = 0. In what follows, the RPA equations for hfe//, Eq. (|^), are obtained. One- and 
two-body sectors of hfe//, Eq. (^), specified in Eqs. (0). o and (^), contribute to 
the RPA. Calculating the commutators of Eq. (|103|) with the effective Hamiltonian H^ff, 
the RPA equations are obtained as 


M„X"'’(fc) = 2e(k)X‘d(k) _ 2 J ^i;d^(k,q)X‘-'(q) - 2 j ^i;^‘-'(k,q)Y‘-'(q) 

- 2 j ^GfJ-’(k,q)X‘-'(,q) - 2 j ’(fc, ,)y‘’(g) 

Af„r""(fc) = -2E(k)Y‘‘^(k) + 2 1q)Y‘''{q) + 2 / ’(*=. <l)X‘-’{q) 

+ ^ + 2 / qXW), dos) 

where an effective single particle energy e{k) is dehned in Eq. (|6^) . Here, the instanta¬ 
neous tensor terms I{k,q) are 

Cx^''ik,q) = CfVL+c{k,q){vl{-k)v^{-q)){ul{q)u4k)) 

= -CfVL+cik,q){vl{-k)usiq))ivl{-q)Uaik)) (106) 

Iyx^^{k,q) = -CfVL+c{k,q){ul{k)v^{-q)){u\{q)vp{-k)) = Ily, 
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and the generated terms G{k, q) are 

q) = CfWi{k, q){vl{-k)aiV^{-q)){u\{q)ajUa{k))D,j{k - q) 
q) = GfW^ik, q){ui{k)aMq))H{-q)aM-kmAk - q) = Gl 
G'^^^^{k,q) = -GfW2{k,q){vl{-k)aiUs{q)){vl{-q)ajUc,{k))Dij{k - q) (107) 
Gy^G(^k, q) = -GfW2ik, q){ui{k)aiV^{-q)){ul{q)ajV/ 3 {-k))Dij{k - q) = G^y , 


with potential functions given by Eqs. (^) and Using Eq. ( 0 . we represent the 
instantaneous terms, Eq. (|106|) , as 


K^‘\k,q) = 


+ 



+ 


lT"\k,q) 

GfVL+c{k, q)^ [(1 + s(fc))(l + s{q))6f}^6sa 

(1 - s(fc))(l - s{q))xl(T 2 (T ■ k(T ■ qa2X^xW ' ' ^Xa 

c{k)c{q){6saxl(^2(T • k(T ■ qa2X'y + ■ kxa) 

GfVL+c{k, q)- [(1 + s{k)){l - s{q))xl(T 2 fT ■ qx5x\(^2fT ■ qXa 
(1 - s(fc))(l + s{q))xlcT 20 - • kxsx\(^ 2 (T ■ kXa 
c(fc)c(q)(x^o-2cr ■ qxsxl(^2(y ■ kXa + xl(^20- ■ kxsxl(^2(y ■ qXa) 


(108) 




and use them in further calculations. Analogous expressions can also be found for the 
generated terms G{k,q), Eq. (|107|) . A crucial test of any approach dealing with chiral 
symmetry is the ability to describe the pseudoscalar meson channel. We therefore consider 
an application of Eq. (|105|) to tt, = 0’*’+ L = S' = 1,J = 0 pseudoscalar, and p, 
= O”'*' L = S' = J = 0 vector, states. Based on the quantum numbers of these 
states, the tensor structure of vr- and p-wave functions can be identihed as 


Xf{k) = {ia2TPX^{k),Yf{k) = {-ia2rPY^{k) 

Xf{k) = {aia2T^Xy{k ), Yf{k) = {-ia2CTr%{k ). (109) 


The normalization is chosen 


/ (V(fc)V(fc) - Y'(k)Y(k)) = 1, 


reducing the normalization costant of the full wave function, Eq. (|102|) , to 


( 110 ) 


A^ = !?/>„) = 277c, (111) 

where the factor 2 comes from the trace in the spinor space, and Nc = 3. The RPA 
equations for the momentum wave function components X{k), Y(k) have the same form 
for TT and p states 


MnX{k) = 2e{k)X{k) - j ^ I^^{k, q)X(q) - I^y{k, q)Y{q) 
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- / ,)y (,) 

-M„F(fc) = Mk)Y(U) - / ,)y(,) - / ,)-y(,) 

- /^^G,,(fc.,)y(,)-/^^G,.(fe.,W,). (112) 

where the kernels I and G for vr are 

Ilx{k,q) = Iyy{k,q) = CfVL+c{k,q)]^[{l + s{k)){l + s{q)) 

+ (1 - s(fc))(l - s(q)) + 2c(fc)c(g)a;] 

Ily{k,q) = Iy^{k,q) =CfVL+c{k,q)]^[-{l +s{k)){l-s{q)) 

- {1 - s{k)){l + s{q))+ 2c{k)c{q)x] 

G:,(k,g) = G-Jk,q) = 2CfW,(k,g)^ - (1 + s(fc))(l - s),)) 

- (l-.(fc))(l + .(,))-2.(fc)c(,)<i±ii!^fi^^‘ 
G',(fc,9) = Gyfc,,)=2C/H'2(A;,,)i (1 + s(/!))(l + s),)) 

+ (1-.(/=))(!-.(,))-2c(/b).(,)h±£!l^^ihl±ih]. (n^ 

and for p are 

IL{k,q) = I^y{k,q) = CfVL+c{k,q)]^ {1 + s{k)){l + s{q)) 


Ixyik,q) 

GL{k,q) 

G^y{k,q) 


^(1 - s(fc))(l - s{q)){4x‘^ - 1) + 2c{k)c{q)x 

q) = C^f4+c(fc, 9)^ [ -^(1 + <k)){l - s{q)) 
i(l - s{k)){l + s{q)) + ‘^c{k)c{q)x 
G^yyik, q) = GfW,{k, q)^ 1(1 + s{k)){l - s{q)) (^1 - 

i(l - s(k))(l + .(,)) (l - - l<k)c(g) (x + 

G;,(fc, g) = C,W,{k, g)l [i(l + s(fc))(l + s{g)) 

1(1 - s(k))(l - s{g))(2x'‘ - 1) + |c(fc)c(^) (x - j 


where we introduced x = k ■ q, and used 




(1 — x^)k^ 

{k - qr)2 


26 



k iq i 


x{k‘^ + g^) — (1 + x‘^)kq 


(115) 


(fe - qf 

with I = k — q. The obtained vr and p RPA eqnations are IR hnite as discnssed in the 
following. Consider the collinear limit k ^ q for the instantaneons terms, I{k, q), since 
the conhning potential in Vl+c canses the IR problem. We get 

VL+cik, 9) [ 1 + s{k)s{q) + c{k)c{q) ] ^ 2(fc - q)-^ 

Kx yL+c{k, q) [ 1 + s{k)s{q) + c{k)c{q) ] ^ 2{k - q)"^ , (116) 

and 


Ily VL+c{k, q) [ -(1 - s{k)s{q)) + c(fc)c(q) ] 

^ VlMK q)0{k - q) + 0{{k - q)-2) 

Ixy VL+c{k, q)^ [ -(1 - s{k)s{q)) + c{k)c{q) ] 

^ VL+c{k, q)0{k -q) + 0{{k- q)-^) . (117) 

The TDA kernel I^x has the same IR behavior for both channels and is IR singnlar. 
However the effective energy e{k) behaves the same way in the IR, Eq. (|80|) , and comes 
with an opposite sign in the TDA/RPA eqnation, cancelling exactly the IR divergence. 
The IR behavior of the TDA kernel shonld be same for other channels, and can be nsed 
to check calcnlations. In the RPA kernel, Ixy, the hrst term 0{k — q) disappears after the 
angnlar integration, and the second term 0{{k — q)“^) converges in the intergral. Thns, 
the TDA and RPA eqnations are IR hnite for the conhning potential. In the UV the 
potential part of interaction, which contains intergrals with kernels I and G, is regnlated 
by the wave fnnctions X and Y, vanishing for large momenta. The kinetic part contains 
the UV hnite ehective energy, e{k), Eq. (|6^, which has been renormalized by adding 
the mass connterterm. 

The RPA eqnations are the eigenvalne problem for which can be diagonalized in 
the twice size space of {X, Y), compare to the TDA reqniring the size of only X. In the 
matrix form, Eq. (|112|) is written as 


with 


A{k,q)X{q) + B{k,q)Y{q) = 
-B{k,q)X{q) - A{k,q)Y{q) = 


MX{k) 

MY{k), (118) 


A{k,q) = {2e{k)6f^ q-Fxx{k,q))dq ^ 2e{k) - dqFx;x{k,q) 

B{k,q) = -dqFx;y{k,q), (119) 


where F inclndes the instantaneons and generated terms, Eqs. ( |113|) and (|114| ); Fxx{k, q) = 
Ixx{k, q) + Gxx{k, q), and the same for xy component. In Eq. (|118|) the integration over 
q is implied, and we omited factors (27r)^. The RPA matrix size can be rednced by a 
factor of 2. Using variables 


= X±y, 


( 120 ) 
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the RPA equations are given by 


{A + B)'il)+ = M%Ij_ 

= Mij+, ( 121 ) 


which can be decoupled, at the expence of more complicated kernel, as 
[{A-B){A + B)]{k,q)Mq) = M^^k) 

[{A + B){A-B)]{k,q)i:_{q) = . (122) 


The diagonalization of either [A — B){A + B) or [A + B){A — B) gives the eigenvalue 
M for the RPA problem. Taking R = 0 we come back to the TDA equation. 

Finally, we calculate pion decay constant in the TDA and RPA schemes. Using Thou- 
less’ theorem applied to the chiral charge in one flavor case 

Qs = J dxi)\x)'y5'i/j{x), (123) 

and the effective Hamiltonian, one gets 

{Q\[Q,,[Q,,H,ff]]\Q) = (124) 

Using I = J2n iV'n) (V'nl/A^, the l.h.s. can be written as 

{n\[Q,,[Q,,H,ff]]\n) = - 2 ^i|(D|g 5 |^n)rM„, (125) 

n 

where = En — Eo, ife//|'0n) = -Unit/'n) and Heff\fl) = Ro|fl)- Here, the normalization 
constant is At = J2n and the volume is set to unity, U = 1. Dehning the weak decay 
constant as 

fn = , ( 126 ) 

we obtain the Gell’Mann-Oakes-Renner relation 

= -2m,onsSi^), (127) 


where rriconst is the quark constituent mass, and (V’'*/’) is the quark condensate. Using Eq. 
(^), the chiral charge is given by 

Qs = Y J+dl{k)ds{k)] 


s{k) [bl{k)dl{-k)+ds{-k)K{k)]) . 


(128) 


Taking the matrix element of the chiral charge, Eq. 
the RPA pion wave function, Eqs. ( |^ and ( |100|) , we have 


, between the vacuum state and 


u = 


2Nr, 


dk 


J (27r)3 


s{k) {X^{k)-Y^{k)) , 


where from Eq. (|111|) N = 2Nc. For Nc = 3, the pion decay constant is given 


u = 


a/6 f k'^dk 




27r2 


s{k) (x^(fc) - y;(fc)), 


(129) 


( 130 ) 


in the RPA. In the TDA, U = 0 and M-j^ is the eigenvalue of the TDA equation. 













4 Numerical results 


In this section we obtain the numerical solutions of the quark gap equation (subsection 


For the numerical calculations we have used the routines from the SLATEC linear 
algebra archive, part of the Netlib database maintained by UTK and ORNL. The routines 
are found at www.netlib.org/slatec/lin/ and a description of the entire SLATEC archive 
can be found at www.netlib.org/slatec/toc. 

4.1 Energy gap, quark propagator and chiral condensate 


LT) and the TDA/RPA bound state equations (subsection B, and discuss the results. 


Mo, 2kas 
(GeV) 



0 0.2 0.4 0.6 0.8 1 

k{GeV) 


Figure 1: The numerical solution of dynamical quark mass, Mo(fc), and the scalar part 
of the propagator, 2kas, in the chiral limit with conhnement. The parameters for the 
numerical solution of the gap equation are a = O.lSGeV'^, A = IGeV. The results are 
compared with the ht function given by Mo(fc) = 0.0024/(A:^[ln(/c^/0.04)]*’'^^) (parameters 
are in powers of GeV). 

We have numerically solved the linearized gap equation, Eq. (|5^) , for 

M{k) = M^^\k) + 5M{k), (131) 

where the mass correction is denoted as 6M{k), while the mass gap M^^\k) is the hrst 
iteration obtained by the gauss algorithm |P], providing an initial guess. The linearized 
gap equation has the form 


J dqA{k,q)6M{q) = B{k), 


(132) 
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k{GeV) 

Figure 2: The numerical solution of dynamical quark mass, MQ{k), and the scalar part 
of the propagator, 2kas, in the chiral limit when Coulomb and generated potentials are 
added with the running coupling as{k‘^) (same parameters as in Fig.l). The results are 
compared with the £t function given by Mo(fc) = 0.0060/(A;^[ln(/c^/0.04)]°'^^). 


and, using Eq. (PD for sine and cosine, it reads 

SM{k) \k+ f ^ ( C!V,,c(k,q)qx 
j 47r2 + 

^ CfW{k,q)q {k^ + q'^)x - kq{l + x^) /^2 

^g2 + (M(0)(q))2 


q^dqdx 


5M{q) 


47r2 
CfWik,q) 


{k - qy 
CfVi+cik, q) 




^q2 + (M(0)(g))2 q^ + {M^^\q)) 

/ [M^^\q)Yk {k)q + q^)x — kq{l + x"^ 


-fc 


{k-q) 


^q2 + (M(0)(q))2 \q^ + {M^^\q)f q^ + {M^^){q)f 

-/t(Af'">(fc) - m(A)) + / fpp. \ P/U+c(fc, - M‘-"\k)qx) 

(M<«>(q)k - Af(»)(fc),<k±khlTih±£!l') e-V/A" 


^q^ + (Mmq)Y 


{k - qy 


where x = k ■ q, the potential functions Vl+c and W are dehned in Eq. (0)> and the 
running mass m(A) is given by Eq. (p|). 


Acd ! Is? 


( 133 ) 
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Figure 3: One particle dispersion relation, Eo(k) = \JkP‘ + Mo(fc)^, free dispersion, 
Eo{k) = fc, and the vector part of the propagator, l/2cr^, in the chiral limit with con- 
hnement and with conhnement plus perturbative potentials (i.e. when Coulomb and 
generated potentials are added). The parameters are same as in Fig. 


As discussed before the obtained gap equation is free from IR and UV divergences, 
and the mass gap, M(fc), exists for a conhning and Coulomb potentials provided both 
instantaneous and generated by flow equations terms are taken into account. 

In numerical calculations of Eq. 


we have used the hie dgeco.f followed by 
dgesl.f from the SLATEC linear algebra archive, which solves A x X = B. 

For a hxed cut-off A = 1 GeV, the mass gap, Moi^k), and the product of momentum 
and scalar part of propagator, /c(j°(fc), in the chiral limit are displayed in Fig.|l| (conhning 
potential) and in Fig.^ (Coulomb and generated terms are added). Hereafter, Mq denotes 
the solution of the gap equation in the chiral limit, not be confused with the hrst iteration 
in our calculations, Since mA ^ 0, the nonzero mass gap Mq = mA/{l + B) is 

generated, and the chiral symmetry is broken dynamically. We denote the maximum value 
of the gap as Mo(0) and call it constituent quark mass. In the chiral limit constituent 
quark mass is generated as Mq{0) ~ 74 MeV with conhning potential, and slightly bigger 
Mo(0) ~ 95 MeV when perturbative terms are included. At high momenta, from Eq. 
(1^, 2ka^g{k) —> Moi^k), and since BCS pairing is a low-momentum ehect, the mass gap 
function Mo(fc) vanishes rapidly as Mo(fc —>■ cxo) ^ l//c^ —> 0. The behavior of the A{k) 
and B{k) amplitudes in the perturbative or UV asymptotic region is well known from 
the QCD renormalization group and operator product expansion and QCD sum rules 
IT[. The behavior has been summarized in the course of explicit numerical solutions 


and model building [T^. For k"^ ^ ^qcd f^e leading-log result for the chiral mass gap. 
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Figure 4: Sine and cosine of the Bogoliubov-Valatin angle in the chiral limit with con- 
hnement and with conhnement and with conhnement plus perturbative potentials (i.e. 
when Coulomb and generated potentials are added). The parameters are same as in Fig. 

0 - 


Mo(fc) = mA{k)/{l + B{k)), is 


Mo(fc 


K 


OO = 


F[ln(fcVA^^^)]i 


—d 


(134) 


where d = 12/(33 — 2Nf) is the anomalous dimension of the mass, Nf is the number of 
quark flavours, Aqcd ~ 0.20 GeV is the scale parameter of QCD, and k is a constant 
given by 


- 3 ^ ^ 

with the scale = 1 GeV^. We htted our numerically obtained solution, Mo(fc), with the 
function Eq. ( |134| ) at high momenta. For Nf = 6, 1 — d = 0.43, we obtain k = 2.4: MeV 
and K = 0.0 MeV, which corresponds to the chiral condensates {'ip'ip)o ~ —{85.3 MeV)^ 
and ('i/'0)o ~ —(115.9MeE)^ for conhning and conhning+perturbative potentials, re¬ 
spectively. These estimates for condensates are lower than calculated directly from Eq. 
m (see below). Though the condensate is a measure of chiral symmetry breaking, it is 
sensitive to the UV region, and therefore it is enhanced by generated terms from the flow 
equation. Note that we reproduce a correct high momentum behavior predicted from the 
perturbative renormalization group analyses. It happens only if the terms from the flow 
equation together with the Coulomb interaction are added (Coulomb interaction alone 
does not work). 
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Figure 5: Scalar, ag, and vector, a^, parts of the propagator, in the chiral limit with 
conhnement and with conhnement and with conhnement plus perturbative potentials 
(i.e. when Coulomb and generated potentials are added). The parameters are same as in 

Fig. 0. 


Similar behavior of the mass gap M{k) and scalar part of propagator is seen for the 
light M, d quarks, with the current mass m = 8 MeV. The predicted leading-log behavior 
for high momenta 


M{k —>■ oo) = 


m 


[l/2\n{k? / Is?QQjy)Y 


(136) 


provides a htting function, where m is the renormalization point independent current 
quark mass. For Nf = 6,d = 0.57, we obtain m ~ 20 MeV and m ~ 32 MeV for conhning 
and conhning-|-perturbative potentials, respectively. Since m 7 ^ 0, there is an explicit 
chiral symmetry breaking. The htting functions, Eqs. ( |134|) and ( |136|) , were obtained 
in Ref. |15] using covariant calculations. In our calculations, we obtain the same htting 


functions. 

The quark energy dispersion, Eo{k) = \Jk? + M^^k), where Mq is the numerical solu¬ 
tion of the gap equation, and the inverse of vector part of propagator, l/ 2 cr°(fc) = Q{k)/{1 + B{k)), 
in the chiral limit at A = 1 GeV are depicted in Fig.0. The free behavior, Eo{k) = k, 
is recovered at high energies while for low energies constituent quark masses at /c = 0 , 
roughly 36.8 MeV and 46.7 MeV are obtained for conhning and conhning-|-perturbative 
potentials, respectively. The inverse of vector part of propagator exactly reproduces the 
energy dispersion as expected from Eq. (O), l/ 2 ff„“(fc)=E„(fc). 

The sine and cosine of the Bogoliubov-Valatin angle in the chiral limit are depicted 
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Figure 6: Cut-off dependence of the constituent quark mass in the chiral limit (same 
parameters as in Fig.l). Crosses represent solution with confinement. Boxes [diamonds] 
represent solution when Coulomb and generated potentials are added with the constant 
value of coupling gg [with the running coupling 


in Fig.^. At low momenta the cosine behaves linearly, with the slope 1/M(0), 

1 


c{k 0) 


M(0) 


k. 


(137) 


The obtained values of Mo(0) in this way are the same as extracted from the dispersion 
relation. From Eq. (EIP’ the relations s{k) = 2as and c{k) = 2kay give additional insight 
into the behavior of scalar and vector parts of propagator. 

The scalar, cr°, and vector, a^, parts of the propagator in the chiral limit are presented 
in Fig.|^. From Eq. (pT]) , at high momenta both amplitudes vanish rapidly. 


(Ts{k 

ay{k 


oo 


oo 


IM(fc) 

2 k 
1 1 

2k' 


while for low momenta both tend to be constants, 

1 


as{k —> 0 ) 
ay{k ^ 0 ) 


2 

1 1 


1 1 

2 F 


2M(0) 


(138) 


(139) 


As expected, adding Coulomb and flow equation terms increases the scalar amplitude and 
decreases vector amplitude at low momenta, that amplihes dynamical chiral symmetry 
breaking with larger Mo(0) and 
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Figure 7: Quark condensate cut-off dependence in the chiral limit (same parameters as in 
Fig.l). Crosses represent solution with confinement. Boxes [diamonds] represent solution 
when Coulomb and generated potentials are added with the constant value of coupling 
Qs [with the running coupling 5's(g^)]. 


The sensitivity of the constituent quark mass (maximum mass gap), Mo(0), and the 
quark condensate in the chiral limit to the cut-off A is displayed in Fig.^ and |^, respec¬ 
tively. Remarkably, after roughly A = 4 GeV the mass gap saturates to a constant value 
even when Coulomb and generated potentials are added (adding only Coulomb causes 
slow logarithmic dependence). This proves that the obtained gap equation is renormal¬ 
ized completely and does not need counterterms in the chiral limit. The constituent 
masses tend to values 63 MeV for confining and 

Mo{0)^70MeV, (140) 


for con£ning-|-perturbative potentials. Taking a running coupling does not change the 
result for M(0). The chiral condensate, calculated using cut-off independent mass gap 
Eq.(^), rises logarithmically with cut-off A when perturbative potentials are 
added. By including the leading-log nonperturbative running coupling constant [^ 
the gap equation (in the perturbative kernel) 


m 


= 


dn 


ln ((,2 + ■ 


(141) 


with scales /i ~ 0.87 GeV (accounts for freezing of since confinement is present) and 
ka ~ 0.16 GeV, found in |]^ by a fit to the non-relativistic heavy-quark lattice data, and 
d = 127r/(33 — 2A^/), we damp the growth of the chiral condensate. Since combination 
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Figure 8: Dependence of the constituent quark mass on the current quark mass, 
(Mm(0) — m)/Mo(0), where Mm{0) and Mo(0) are the constituent quark masses for the 
current quark mass m and in the chiral limit m —0, respectively (same parameters as 
in Fig.l). Crosses represent solution with conhnement. Boxes represent solution when 
Coulomb and generated potentials are added. 

m{A){'ipip)(A) is renormalization group invariant (it appears in Hamiltonian), the leading- 
log bahavior is given from Eq. (|40|) 

m(A) = m(Ao) - -^^blnA/Ao^ 

W)(A) = (^^)(Ao) - -^^einA/Ao^ , (142) 

where Aq is the renormalization point. Hence the cut-off dependent correction to conden¬ 
sate In A can be absorbed by introducing running coupling g'^ —>■ g‘^{A) ~ 1/lnA. 
Relations between nonperturbative scales Aq, g, and can be found. This procedure 
renormalizes quark condensate to leading-log order and freezes its value at 

(7/;V')o ~-(155Me■F)^ (143) 

while only with conhnement condensate the value of —(92 MeVY (roughly —(100 MeV)^ 
reported in ||^, |^) is obtained. Flow equations improve the chiral condensate by ~ 68%, 
although the obtained value is still low (which is a common feature for most Hamiltonian 
methods). 

Constituent masses and quark condensates for different current masses are shown 
in Fig.^ and |], respectively. In light quark sector (u, d quarks) there is practically no 
deviation from the chiral result. However, constituent masses and condensates increase 
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Figure 9: Dependence of the quark condensate on the current quark mass, 
, where {'tp'ip)m and {'ip'ip)o are the quark condensates for the current 
quark mass m and in the chiral limit m —> 0, respectively (same parameters as in Fig.l). 
Crosses represent solution with conhnement. Boxes represent solution when Coulomb 
and generated potentials are added. 


slowly when approaching the strange quark mass, and subsequently decrease. Generally, 
this behavior is very slow, and depends on particular model used for calculations. In 
the instanton liquid model, it has been reported [|^ that the constitituent mass and 
the condensate decrease as the current quark mass increases. We recover the relation 
between quark and gluon condensates |T^ 


— 1 ry 

m(-V’V’) = , 

K7T 


(144) 


in heavy quark sector. Here the expectation value with respect to the quark vacuum 
state is implied in f{Q\'ip'ijj\Q)f, and the average over the gluon vacuum conhgurations 
stands in g{Q\Ffj,i^F^'^\Q)g, assuming that the quark and gluon vacuum states factorize 
as the direct product |D) = |D)/ ® \^)g- We obtain {^Fg^F^’^) ^ 0.010 GeV^ for strange 


quark, with k = 13.2 |]T^, in agreement with the QCD sum rules 0.012 [0. In the 

light quark sector we obtain too low value for the gluon condensate, suggesting that this 
relation holds only for heavy quarks. 
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4.2 Pion and p-meson bound states 


Numerical solutions of Eq. ( 122 ) in RPA and TDA for pion (Eq. (p.l3| )) and p meson (Eq. 
(|114|)) interaction kernels are obtained variationally with a set of Gaussian test functions, 


i’i 


{n\i) 


n=l,Nt^ 

-k'^nl20^ 


(145) 


where %[)i is an eigenfunction in the momentum space k, discretized by i = 1 ,A^, and 
/3 is a variational parameter, chosen to give the minimum energy expectation value. In 
this way we avoid the direct diagonalization of Hamiltonian in the momentum space, |i), 
where the interaction kernel is infrared divergent, such as Coulomb ~ and conhning 
~ potentials at fc ~ 0. Instead we calculate Hamiltonian matrix elements between 
Gaussian functions, and obtain a regular Hamiltonian matrix in the \n) space, which is 
solved as an eigenproblem. Typically one choses Nh N io achieve convergence and 
stability of a result (large N for numerical integration of interaction kernels Eq. ([I13| ) 
and Eq. (|114D , and small Nh to give several lowest eigenvalues from discrete spectrum). 
Such procedure is important when calculating pion mass spectrum in the chiral limit, 
with the current quark mass m = 0 (since there is no regulator in pion denominator as 
k 0). 



TDA, (MeV) 

RPA, (MeV) 

conf. 

504 

1364 

2115 

222 

1416 

2298 

conf.+Coul. 

608 

1514 

2249 

427 

1521 

2309 

conf.+Coul. +gen. 

513 

1411 

2161 

180 

1413 

2218 


Table 1: Pion spectrum for the ground, hrst and second exited states in the TDA and 
RPA approaches with conhning, con£ning+Coulomb and conhning+Coulomb+generated 
potentials taken. Chiral limit m = 0 (a^ = 0.4, a = O.lSGeR^, A = lOGeV). 


m, (MeV) 

TDA, (MeV) 

RPA, (MeV) 

150 

1038 

1926 

2936 

1037 

1986 

3077 

100 

885 

1762 

2697 

868 

1811 

2826 

50 

716 

1590 

2431 

660 

1626 

2537 

10 

553 

1446 

2212 

366 

1460 

2283 

5 

532 

1428 

2186 

293 

1437 

2250 

0 

513 

1411 

2161 

180 

1413 

2218 


Table 2: Pion spectrum for the ground, hrst and second exited states in the TDA and RPA 
approaches for diherent current masses of constituents. Conhning+Coulomb+generated 
potentials are taken (the same parameters as in the table 1). 


RPA/TDA Eq. (|122| ) has the form MmnXn = XmNmnXn in the \n) space, where 
Nmn = {m\n), Xn and A„ are eigenfunctions and eigenvalues. Note that the RPA matrix 
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TDA, (MeV) 

RPA, (MeV) 

conf. 

659 

1484 

2258 

642 

1482 

2256 

conf.-|-Coul. 

750 

1678 

2515 

732 

1676 

2514 

conf. -|- Coul. -|-gen. 

718 

1592 

2377 

700 

1590 

2376 


Table 3: Spectrum of the p meson for the ground, hrst and second exited 
states in the TDA and RPA approaches with confining, confining+Coulomb and 
conhning+Coulomb+generated potentials taken. Chiral limit m = 0, the same parame¬ 
ters as in the table 1. 


m, (MeV) 

TDA, (MeV) 

RPA, (MeV) 

150 

1130 

2086 

3247 

1128 

2086 

3247 

100 

986 

1916 

2990 

983 

1915 

2990 

50 

839 

1744 

2692 

833 

1744 

2692 

10 

727 

1616 

2436 

714 

1615 

2435 

5 

719 

1603 

2406 

704 

1601 

2405 

0 

718 

1592 

2377 

700 

1590 

2376 


Table 4: Spectrum of the p meson for the ground, hrst and second exited states 
in the TDA and RPA approaches for different current masses of constituents. 
Con£ning-|-Coulomb-|-generated potentials are taken (the same parameters as in the 
table 1). 


Mmn is not hermitian, A x B ^ B x A (i.e. it is not a Hamiltonian) and one cannot 
use algorithms for diagonalization of symmetric matrices. In numerical calculations of 
Eq. ( |122|) we first used the hies dgeco.f and dgesl.f to hnd the product N~^ x M, 
that reduces eigenvalue problem to N~^MX = AX, and then used the hie rg.f from 
the SLATEC archive to solve the general matrix for eigenstates. Results for the pion 
and p-meson masses in RPA and TDA approaches are presented in Tables In the 
chiral limit RPA gives ground state pion masses which are signihcantly lower than those 
obtained using TDA. Including the generated interaction terms reduces ground state pion 
mass even more, increasing the mass splitting between the pion and p meson (Tables ^ 
^). In the chiral limit we get 


= 180 MeV , Mp = 700 MeV 

Mp-M^ = 520 MeV . (146) 


Although Coulomb and generated interactions are both perturbative, they act in dif¬ 
ferent directions: Coulomb [generated] term increases [decreases] a meson mass. Indeed, 
this is in an accordance with the standard perturbation theory where the leading order 
perturbative corrections always shift the energy of a ground state down. Ehects of using 
the RPA instead of the TDA including the generated terms are not so pronounced for 
excited states. 

The TT — p mass splitting of 520 MeV in the chiral limit is close enough to the lattice 
data splitting of 600 MeR [^]. However, we are unable to get zero mass pion either in 
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the BCS or adding the corrections from the leading order flow equations. Various reasons 
why the zero mass pion solution cannot appear within the BCS approach are summarized 
in Appendix The underlying reason of failing to produce the Goldstone boson seems 
to be associated with a breakdown of covariance in the BCS model. Including higher 
orders of calculations is necessary in order to approach to a covariant result. 

In the chiral limit the tt— p mass splitting, = Vfp—with all terms contributed 

is 5Et^p = 205 MeV in the TDA and = 520 MeV in the RPA. These values should be 
compared to SE^^p = 155 MeV in the TDA and hi?,rp = 420 MeV in the RPA, calculated 
with static conhning potential alone. Flow equations improve the n — p mass splitting 
by 32% in the TDA and by 24% in the RPA. 



Figure 10: n — p mass splitting for the ground state in TDA and RPA approaches. 
Chiral limit, m = 0, and confining+Coulomb+generated interactions are taken (same 
parameters as in table[^). Fit function is 12000/(m+20). 

For nonzero current quark masses, the main contribution to meson masses is coming 
from the solutions of the gap equation, while the contributions from bare mass kinetic 
terms are negligible (Tables §). This supports an idea of the constituent quark model, 
where interactions give rise to massive constituents which build a meson mass in a 
valence approximation. Since a numerical solution of the gap equation carries its mass 
dependence, it is difficult to analyze analytically the dependence of the meson spectrum, 
in particular of the tt — p mass splitting, on the current mass. 

Numerically obtained dependence of the tt — p mass splitting, Mp — M^, as a function 
of the bare mass of one of the quarks is shown in Figure |^. We hnd 1/mconst behavior, 
where rriconst is a constituent quark mass, which is valid for heavy quarks and continues 
to be valid for lighter constituent quarks. This fall-off is more rapid in the RPA than in 
the TDA. From the RPA £t function, a constituent quark mass can be approximated. 
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1.09 


M(A)/M(10) 



Figure 11: Cut-off dependence of the pion mass in the chiral limit M(A), normalized to 
the pion mass at A = lOGeV^, M(10), (same parameters as in table|^. Line with crosses 
represents RPA solution with con£nement-|-Coulomb, the line with boxes stays for RPA 
solution when generated potentials are added. 


uniformly for heavy and light quarks, as rriconst = m + 20 [MeV), with the bare quark 
mass m. The 1/rriconst behavior is characteristic for the hyperhne interaction. However, we 
reproduce this behavior taking into account both spin dependent interactions and vacuum 
effects of the dynamical chiral symmetry breaking. We do not separate the effects of the 
hyperhne interactions and the chiral symmetry breaking, since after Bogoluibov-Valatin 
transformation all quark interactions are calculated in a chiral noninvariant framework. 

We reproduce 1/rriconst behavior throughout the whole range of current quark masses. 
One may conjecture that the hyperhne interaction dominates for heavy quarks, and the 
conhning potential and the chiral symmetry breaking are dominant for light quarks. 
However, we do not see a change in regime. We thus conclude that I/rriconst behavior of 
the 71 — p mass splitting, attributed to the hyperhne interaction, is actually due to both 
hyperhne and the chiral symmetry breaking ehects, which is consistent with the lattice 
calculations and experiment |]^ 


Formally, the static interactions of the form which include the sum 

of conhning and Coulomb potentials and initiate the breaking of the chiral symmetry, is 
referred as the chiral symmetry breaking interactions, and the generated terms, having 
the structure are referred as the hyperhne interactions. 

We hnd in the chiral limit that in the TDA roughly 30% of the tt — p mass splitting 
is due to the presence of the hyperhne interaction and the rest 70% is due to the chiral 
symmetry breaking. In the RPA this ratio is 40% for the hyperhne and 60% for the chiral 
symmetry breaking. However, the numerical value of this ratio depends on the details of 
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the confining interaction. In our calculations we used the linear rising potential generally 
accepted in the quark model phenomenology with the string tension a = O.lSGeld^ 
predicted by the lattice studies. It is important to have both terms to reproduce a 
correct mass splitting. 



fc(GeV) 


Figure 12: Ground state pion wave function in TDA and RPA approaches. Chiral limit, 
m = 0, and confining+Coulomb+generated interactions are taken = 0.4, a = 
O.lSGel/^A = WGeV). 

The dependence of the pion mass in the chiral limit on the cut-off parameter, M(A), 
is shown in Fig.[TI|. The RPA solution, with only confining-l-Coulomb potential included, 
grows unlimited (due to Coulomb), while adding the generated term stabilizes M(A), 
which saturates roughly at M(10). Stable result confirms that the TDA/RPA equations 
are completely renormalized when the generated terms by flow equations are included. 

The pion wave functions for the ground and first excited states are depicted in Figs.|T^ 

between the TDA wave function, ipTOA, and the creation component of the RPA wave 
function, X^pa, is most significant for the pion ground state, while it is smaller for the 
pion excited states and for the p-meson ground state. We conclude, as expected, that the 
pion ground state is most sensitive to the improvement from the RPA over the TDA. 


13. The ground state p-meson wave function is also shown in Fig.|l^. The difference 



TDA, (MeV) 

RPA, (MeV) 

conf. 

20 

45 

conf. -I-Coul. -|-gen. 

39 

92 


Table 5: Pion decay constants in the TDA and RPA approaches with confining and 
confining-|-Coulomb-|-generated potentials taken. 
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Figure 13: Pion wave function for the first exited state in TDA and RPA approaches. 
Chiral limit, m = 0, and con£ning-|-Coulomb+generated interactions are taken into 
account (same parameters as in Fig. [T0|) . 



fc(GeV) 


Figure 14: Ground state p meson wave function in TDA and RPA approaches. Chiral 
limit, m = 0, and conhning+Coulomb+generated interactions are taken into account 
(same parameters as in Fig.pRf). 
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Using the TDA and RPA wave fnnctions, we obtain the pion decay constants given 
in the table We find higher valnes of in the RPA than in the TDA. Taking the 
confining potential we obtain = AbMeV in the RPA, that is higher than given in 
Ref. 0 /tt = 20MeU. Adding the dynamical terms improves the pion decay constant, 
/tt = 92MeV, that should be compared with the experimental value /jr = 96MeV. The 
Gell’Mann-Oakes-Renner relation, Eq. (|127|) , is satisfied within 6%, showing how close 
our wave functions are to the exact Hamiltonian eigenfunctions. This may explain why 
we obtain a realistic pion decay constant, 92MeV, in the RPA, although we obtained 
low value of quark condensate —(155MeU)^ in subsection 


5 Summary and Conclusions 

Starting with the leading order QCD Hamiltonian in the Coulomb gauge and introducing 
phenomenological potential to account for the nonperturbative effects, we performed a 
sequence of transformations of the Hamiltonian: BV transformation to a massive quasi¬ 
particle basis and integral transformation achieved by flow equations. After BV transfor¬ 
mation the phenomenological potential becomes strong, saturating strong QCD interac¬ 
tions, and describes the BCS solution, which involves dynamical breaking of the chiral 
symmetry. Residual interactions (weak in the quasiparticle basis) are treated perturba- 
tively by flow equations and provide perturbative corrections to the BCS solution. In this 
way we have utilized the scheme of the BCS model, where nonperturbative features such 
as dynamical chiral symmetry breaking and massive quasiparticle modes are explicitly 
present, and included perturbatively dynamical interactions in the BCS framework. 

Including dynamical interactions by flow equations, we find the correct ultraviolet 
behavior of Hamiltonian solutions. Namely, for the first time the renormalized gap and 
Bethe-Salpeter equations which are finite in both UV and IR regions are obtained. More¬ 
over, no additional UV renormalization is required in the chiral limit. 

Adding dynamical interactions only slightly enlarges the quark dynamical mass; 
mo(0) = 70MeV in the chiral limit. However, the chiral condensate has been improved 
by 68% using flow equations and equals ('0V’)o = —(155 MeU)^. This value is still low 
compared to the predicted lattice value = —(250 MeV)^, but such feature of lower 

value is common for all Hamiltonain methods. 

The cut-off dependence of the quark condensate at nonzero current masses has been 
obtained which provides the product m(A){'0'^)m(A) invariant with respect to A. The 
chiral condensate {'4>'ip)o is renormalized by including the renormalization group running 
of the strong coupling constant as{q^)- 

Bethe-Salpeter equation is solved in the TDA and RPA, which approximate the 
pion as a valence qq pair in the TDA and allowing in addition to a qq creation for a 
qq annihilation operator in the RPA. The pion ground state mass is most sensitive to 
the improvement over the TDA approach from the RPA. As expected, the pion ground 
state is obtained lower in the RPA than in the TDA. However, we are unale to get 
zero mass pion either in the BCS or adding the leading order corrections from flow 
equations. An underlying reason of failing to produce the Coldstone boson might be 
a breakdown of covariance in the BCS model, which can be cured by including higher 
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orders of calculations. Indeed, even in including the leading order dynamical interactions 
in the BCS solution shifts the RPA pion mass of the ground state down, as expected 
from the standard perturbation theory. We obtain for the RPA pion and p-meson ground 
states = 180 MeR and Mp = 700 MeR, respectively, in the chiral limit. 

In the chiral limit the tt — p mass splitting with all terms contributed is 5M.^p = 
205 MeR in the TDA and hAArp = 520 MeR in the RPA. These values should be com¬ 
pared to SMt^p = 155 MeV in the TDA and = 420 MeR in the RPA, calculated 
with the static conhning potential alone. Flow equations improve the n — p mass splitting 
by 32% in the TDA and by 24% in the RPA. 

One may conjecture that the hyperhne interaction dominates for heavy quark masses, 
and the dynamical chiral symmetry breaking is dominant for light quarks. However, we 
cannot separate these two effects, distinguishing different regimes. By switching one 
contribution or another, we hnd in the chiral limit, that roughly 30% of the n — p mass 
splitting in the TDA is due to the presence of the hyperhne interaction and the rest 70% 
is due to the chiral symmetry breaking. In the RPA, this ratio is 40% for the hyperhne 
and 60% for the chiral symmetry breaking. However, the numerical value of this ratio 
depends on the details of the conhning interaction. In our calculations we have used the 
linear rising potential generally accepted in the quark model phenomenology with the 
string tension a = 0.18 GeR^ predicted by the lattice studies. 

It is crucial that the how equations generate the dynamical spin-dependent hyper¬ 
hne interactions which depend on all in- and out-going quark momenta as well as the 
momentum of propagating dynamical gluon and include the momentum dependent so¬ 
lutions of the gap equation. In the NJL type model as shown in the second paper of 
Ref. 0, the static hyperhne interaction contribute a dominant part to the tt — p mass 
splitting, indicating that this model is close to the nonrelativistic CQM. However, our 
model indicates that dynamical interactions in the BCS framework need to be included 
using how equations to improve the CQM. Thus, the how equations seem to provide a 
working tool to incorporate the corrections toward a covariant result. 


Acknowledgments. The author (E. G.) would like to thank Chueng-Ryong Ji for 
reading the manuscript and providing grammatical corrections. 


A Complete QCD motivated Hamiltonian 


Here we specify, in second quantized form, the full QCD motivated Hamiltonian in the 
Coulomb gauge expanded in the basis of Eq.(|]). We ignore the pure gluon non-abelian 
terms, since they do not contribute to the quark sector. One-body operators and conden¬ 
sate terms arise from normal ordering with respect to the trial vacuum state ID). The 
upper index over a Hamiltonian operator (e.g. denotes the order (power) in 

coupling constant. 

We specify the y-matrices, that are used in this work. 


li 


f 0 ai 

V -Ui 0 


/? = 7o 


1 0 
0 -1 


, 0^2 (^'li 


0 

CTi 

\ 

/ 0 

1 A 

CTi 

0 . 

1 , 75 = 

1 

0 ) 


45 


( 147 ) 


and the spinors are defined in Eq.(^). 

Free quark and gluon part: Eq.(|) inclndes the qnark Kg and glnon Kg kinetic energies, 
K = Kq + Kg, and the corresponding condensate terms, O = Og + Og. The qnark kinetic 
term reads 


(ul{k){ak + f5m)us{k)bl{k)bs{k) — vl{—k){a.k + f3m)vs{—k)dl{—k)ds{—k) 

+ ul{k){ak + Pm)vs{—k)bl{k)dl{—k) + vl{—k){a.k + (3m)us{k)ds{—k)bs{k)^ , 


which, using the spinors in Eq. m and 7 -matrices in Eq. ( |147| ), is reduced to 

dk 


K„ = 


H J {{kc^{k) + rns^{k))[bl{k)bs{k) + dl{k)d,{k)] 

+ (fcs$(fc) - mc$(fc))[6j(fc)dj(-k) ds{-'k)bs{k)]j 

= 11 / + ^2 (ce{k)[bl{k)bs{k) + dl{k)ds{k)] 

+ se{k)[bl{k)dl{-k) + ds{-k)bs{k)]j , 


(149) 


where m is the bare quark mass, and we have used the connection between the nonper- 
turbative angle $ and the Bogoliubov-Valatin angle 0 (See Eq. (|205|) in Appendix 0). 
In order to separate the zeroth and second order contributions we represent the quark 
kinetic term, using Eq. (0) , as 


Kg = 0 / + H{k)bs{k) + dl{k)ds{k)] 

+ 0 / [s^{k)[bl{k)bs{k) + dl{k)ds{k)] 

- c^{k)[bl{k)dl{-k) + ds{-k)bs{k)]^ 

= (150) 

where Kjj^^ and K'^^ correspond to the zeroth and second order quark kinetic energy, 
respectively. The gluon kinetic energy is given by 


K, = ^r=E 


dk 

(27r)' 


■uj{k)a’l\k)a‘l{k ), 


(151) 


Summing the zeroth order kinetic energy yields 

= Y.j + + E / ik)al{k), 

(152) 
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with the effective quark and free gluon energies, E{k) = + M'^{k), Eq. ( 0 . and 

uj{k) = /c, respectively. The quark condensate reads 


= KV 


/ rj JU 


(27r 

with the volume V = (27r)^(5*^^^(0). It is reduced to 

dk 


(153) 


O, = -AN,V 


= -ANrV 


(2vr)= 

dk 

(2vr)' 


■ (fcc$(fc) + msi^,{k)) 
■\Jk‘^ + m2C0(fc), 


(154) 


where we have used Eq. (|205|) . Separating the zeroth and second order, one has 


0„ = -ANrV 


dk 


" (27r)3 

= Of+of. 

The gluon condensate is given by 


+ M2(fc) -AN^V 


dk 

W)'- 


(m — M{k))s^{k) 


0, = 0f = («■=-l)v/^c.(fc). 


(27 

Therefore, the zeroth order condensate term is given by 

dk 


( 166 ) 


( 166 ) 


O'”’ = + 


dk 


(157) 


where the volume is V = (27r)3(5(3)(0), and E{k) = \Jk'^ + M^{k). 

Instantaneous interaction: Eq. (^ includes the linear confining and Coulomb inter¬ 
actions in the quark sector 


v;„„ = E / f n 

Si...54 


(27r)^5® (fci + k‘i-k2- fc4)I4+c(fci, ^ 2 ) 


=1 (27r)3/ 

X : [<(fc4)6f fci) +<(-fc0ds,(-fcl)]T“K(fc2)&s.(fc2) +^^s.(-fc2)4(-fc2)] 

[“L (^ 3)^13 (fcs) + v\^{-k^)ds^{-k^)]T^[us^{k^)bsSki) + VsS-ki)d\^{-ki)\. 

(158) 


= _L 

^ inst ' ^ tnst '} 


with Vi+cik, q) —»• VL+c(fc — q)- Here the linear and Coulomb terms are represented by 
and respectively, and VL+c{k) is defined by 


cy CT 

CfVL+c{k) = 2vrOj^ + 4vr^, 


(159) 
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with the fundamental Casimir operator Cf = T“T“ = {N^ — l)/2Nc = 4/3. Terms arising 
from normal ordering are the one-body operator (self-energy) and the condensate term. 
Self-energy is written as 

^inst = n/ (160) 

|^6j(fc)6,(A;) ^[(Mj(fc)M,,(qr))(M^,(qr)M,(fc)) - (Mj(fc)u,/(-q))(u],(-qr)M,(fc))] 

- dl{-k)ds{-k) Y,[Hi-k)usiq)){ul,{q)Vsi-k)) - (ul(-fc)u,/(-q))(4(-qr)u,(-fc))] 

s' 

+ ds{-k)bs{k) Y,[{vl{-k)us>{q)){ul,{q)us{k)) - (uJ(-fc)u,/(-q))(4(-qr)M,(fc))] 
s' 

+ bl{k)dl{-k)Y,[{ul{k)Us'{q)){ul,{q)Vs{-k)) - {ul{k)vs'{-q)){vl,{-q)vs{-k))]\ , 

s' / 

which, using Eq. (p!oD, is reduced to 

^inst = J-I^Cf^L+c{k,q)[bl{k)bs{k)+dl{-k)ds{-k)] 

\ ((1 + s{k)s{q) + c{k)c{q)k ■ q) - {1 - s{k)s{q) - c{k)c{q)k ■ q)) 

+ E/ ^C^/V,+c(fc,g)[&I(fc)4(-k) +4(-k)6,(fe)] 

^ (2{-c{k)s{q) + s{k)c{q)k ■ g)) 

= E / {sik)s{q) + c{k)c{q)k ■ q) CfVL+c{k, q) 

X [b\{k)bs{k) + dl{-k)ds{-k)] 

+ E / (-c(fc)s(q) -1- s{k)c{q)k ■ q) CfVL+c{k, q) 

X [bl{k)dl{-k) + ds{-k)bs{k)]. (161) 

The correction to quark condensate is given by 

Oinst = Nc [ ^I^J2[('^l(-^)^s'{q)){ui{q)Vs{-k))]CfVL+c{Kq). (162) 

•' s.s' 


which simplihes to 

Oinst = 2iVc J (l - s{k)s{q) - c{k)c{q)k ■ q) C'/ 14 +c(fc, q) . (163) 

Dynamical interaction: Eq. (^ includes to order 0{g) the quark-gluon coupling 


- -OY. 
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dkr^ 


S\^S2-,Ci ' 


(2vr)= 


i2Ti)H^^\kr-k2-ki,) 


V^2a;(fc3) 


X : [u^fci)6jj(fci) +vl^{-ki)ds^{-ki)]T^ai[us^{k2)b,^{k2) + Vs^{-k2)d\^{-k2)] 

Kiks)+af{-k3)]: , (164) 
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where the colon denotes normal ordered products, i.e. all creation operators are on the 
left and annihilation on the right. One-body and condensate terms diverge in the UV 
region, since they contain held operator products at one point. They are regulated as 
discussed in Section 


B Second order flow equations 

Solving the second oder how equation for the particle number conserving part Eq. 

we generate three types of Hamiltonian operators: two-body ehective quark inter¬ 
actions, one-body self-energy terms and condensates. We consider each separately. In 
what follows, s$(fc) = s{k) and the same abbreviation is dehned for a cosine. 


B.l Effective quark interaction 


We calculate an ehective quark interaction only in the sectors which contribute to the 
bound state equations: TDA and RPA. Using expressions for the ehective coupling con¬ 
stant, Eq. ( 0 ). and the generator, Eq. (|l^, the following ehective interaction is gener¬ 
ated, which contribute to TDA equation. 


dVgenif) 


dl 


X *2, fci - k2)gi(ki, kt, ki - hi) + ihiki, ki, fc, - ki)gi(ki,ki, ki - * 2 )] 

X :bHki)T%„(ki)dJ-ki)rdH-k,): 

+ Y. fill + ki-ki- k,) 

SI...Si'' \n=l / 

X [ 771 (^ 2 , fci ,-{ki - k2))gi{k4, fca, -(fci - ^ 2 )) 

+ Vi{ki, fcs, -(fci - k2))gi{k2, ki, -(fci - ^ 2 ))] 

X ■.dsA-k,)nl{-k2)bl{ks)T%si{k4y. , (165) 


Integrating this how equation, we obtain the generated interaction Vgeni} —^ cxo) = Vgen 
(the initial value is Vgen{l = 0) =0); 


^gen 


X 

X 


E n 


dkr 


S1...S4 


=1 (27r)= 


(27r)35(3)(fc^ + fc3_fc2-fc4) 


^(kl - kP " Df + D? (ki)ain.,(k.i))(vl,(-ki)aiV.,(-ki)) 
:bUki)T‘‘b,,{ki)d,,{-ki)T‘‘dU-kiy. 
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i: n 


S1...54 ' 


dkr, 




{27r)H^^\ki + ks - k2 - k^) 


X 


Dij{ki — ^ 2 ) 2 -^1 “I" 


2a;(fci-fc2) Df + Df 

X :4i(-fci)T“<(-fc2)6l3(fc3)T“6,,(fc4): . 


vl,{-ki)aiVs2{-k2)){ul^{k3)ajUs^{k^)) 


( 166 ) 


We combine both terms, using symmetry property of the polarization sum Dji{q) = 
Dij{q), and obtain 


K 


gen 


= E 


Hh 


51 ...54 


“1 (2vr)= 


(27r)^(5®(fci + fcs - ^2 - ^ 4 ) 


X 


Dijiki — ^2) 2 -Di + -D'l -Di + -Di ^ ^ I ^ t ^ ^ t 1 \\ 

2u(k, - k,) 0 (sfTw+srr5fj «(fc.)«.«..(fe))«(-fc3)a,«..(-*.)) 


X :6l.(fe,)r“6„(fe2)rf„(-fe3)r“<(-fe4): , 

where the energy denominators are 

Di = E{ki)-E{k2)-uj{ki-k2),D[ = E{k:,)-E{k^j-u{ki-k2) 
Di = E{k2) - E{ki) - uj{ki - k2), D[ = E{ki) - E{k^) - uj{ki - k2) 

An effective interaction, which contribute to RPA, is dehned by 


(167) 


( 168 ) 


dVgeniP) 

dl 


= E 

51 ...54 ' 

Diji^ki 


dk 


X 


(27r)3 
k2^ 


(27r)3(5®(fci + ^3 - ^2 - M 


2ct;(A:i — ^2 


ul^{ki)aiVs2{-k2)){ul^{k3)ajVs^{-k4^)) 


X [?7o(A: 4, fcg, ki - k2)gv{ki - ^ 2 , ^ 2 , ki) - r]r{ki - ^ 2 , ^ 2 , ^ 1 ) 510 (^ 4 , fcg, ki - ^ 2 )] 
X :6^fci)T“4(-fc2)6l3(fc3)T“<(-fc4): 


i: n 


dkr 


(27r)^(5®(fci + ^3 - ^2 - k^) 


S1...S4 • 


X 


Dij{ki - k2 


2u{ki — k 


L\ (2vr)^ 

vl^{-ki)aiUs2{k2)){vl (-k3)ajUs4ik4 )) 


X [?7o(A:i, k2, ki - k2)gv{ki - ^ 2 , k^, k^) - gv{ki - ^ 2 , ^ 3 , ^ 4 ) 510 (^ 1 , ^ 2 , ki - ^ 2 )] 

X ■.ds4{-ki)T%s,{k2)ds,{-k3)T%s4{k,y , (169) 


that gives after integration 


K 


gen 


= 2: n 


dkr 


X 


(27r)3^ 
9^ 


Dij{ki k 2 ) 2EQ ( t 


2u;(fci-fc2)'' Dl + Dl 

X : hi (fci)T“<(-fc2)6l3 (fc3)T“<(-fc4) 


{27rl6^^\ki + ks - k2 - ^ 4 ) 

ul{ki)aiVs^i-k2)){ul{k3)ajVs4i-k4)) 
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( 170 ) 



(27r)V3)(fci + fc3-fc2 


K) 


X 

X 


Dij(^ki k^) 2^0 -^i' 

-9 


2uj{ki-k2)^ D^ + Dl 
■.dsA-ki)T%sAk2)ds,{-k^)T%sAkA)-. , 


vl^ {-ki)aiUs^{k2)) {vl^{-k3)ajUs^ (k^)) 


where the energy denominators are 


Do = -{Ei^ks) + E^ki) + uj{ki - k 2 )), Dr = uj{ki - k 2 ) - E{ki) - E{k 2 ) 
Do = — [E[ki) + 77 (^ 2 ) + uj[ki — ^ 2 ))) Dr = u}(ki — ^ 2 ) — E[k^) — E{k^ 


Here the energy denominators carry the same lower indeces Di as the generators r]i and 
coupling constants gi which they correspond to. One can further simplify the generated 
interactions Eqs. (|167| ) and (p.7(]| ). In the c.m. frame the generated interactions contribut¬ 
ing to TDA {X component of the RPA wave function), Eq. (|167D , are 


= T, f^2Wi{k,q)D„{k-g) (172) 

X (^{ul{q)aiUa{k)){vl{-k)ajVj{-q)): 6j(qr)T“6„(fc)d^(-fc)T“d;^(-q): 

+ iul{k)aiU5{q)){v\{-q)ajVg{-k)): bl{k)T%s{q)d^{-q)T''dl{-k): ) , 

and to RPA {Y component of the RPA wave function), Eq. dno) , are 


(173) 

0/357 ' 7 

X {^{v\{-q)aiUa{k)){vl{-k)ajUs{q)y. d^{-q)T%a{k)dg{-k)T%s{q)\ 

+ {u\{q)aiVfi{-k)){ui{k)ajV^{-q))-. 6l(qr)T“4(-fc)6),(fc)T“4(-qr): ) , 

where potential functions are given by 


CfW^{k,q) 

CfW2{k,q) 


1 __ 

2uj^{k~q) + {E{k)-E{q)y 

1__ 

2a;2(fc-q) + (E(fc) + E(q))2 ' 


(174) 


B.2 Self-energies 


Two contraction terms from the commutator [r]^^\V^p] in Eq. (|2^) contribute to the 
self-energy operators. In the diagonal one-quark sector [b^b and d^d) the flow equations 
are written as follows 


dSgen(0 _ f dkdq^ Dij{k - q) 


= E 


:Cf 


dl J (27r)® ^ 2uj{k — q) 

X ^iu\{k)(^iUs'{q)){ul,{q)ajUs{k))2rji{k, q,k- q)gi{k, q,k-q) 


(175) 


J2{ul{k)aiVs'{-q)){vl,{-q)ajUs{k))2go{k, q,k- q)go{k, q,k-q) \ b\{k)bs{k ), 


51 


















( 176 ) 


since rjQ^p, k, q) = rj^ik^p, q) (the same holds for Qq) and 
d'^geni]^) 


S/ 


dkdq^ Dij{k — q) 


rCf 


dl J (27r)® ^ 2uj{k — q) 

X ( Y.iyl{-k)aiVs\-q)){y\,{-q)ajVs{-k))2qi{k, q,k- q)gi{k, q,k-q) 

- q,k- q)go{k, q,k- q)^ dl(-k)d,(-k) 

Integrating these flow equations with the generators and coupling constants given by 
Eq. (^) and Eq. ([2l|), respectively, produces the second order self-energy correction, 

6E = E(l) — E(lo = 0). The self-energy operator, at the scale A with / = 1/A, is given by 


^ f dkdq^ ^ 2 dd>ij{k q) 

S.en(A) - 2u:{k-q) 


(177) 


X |^^(Mi(fc)aiM,/(qr))(Mj,(qr)ajMs(fc))—e^ 2Dyx-^) 

- Y.^ul{k)aiVs'{-q)){vl,{-q)ajUs{k))^e^-‘^^^o/^‘'^\ bl{k)bs{k) , 

-^0 / 


and 


_ ^ fdkdq 2 Dij{k q 

S,™(A) - A- j Wf - o 


2uj{k — q) 


(178) 


X 


Y,{vl{-k)aiVs'{-q)){vl,{-q)ajVs{-k))^e^ 

K af 1 


- Y^{vl{-k)aiUs'{q)){ul,{q)ajVs{-k))^e^ d\{-k)ds{-k), 

s' ^0 / 


with 


Do = -{E{k) + E{q)+u;{k-q)) 

Dl = E{k) - E{q) - uj{k - q) . (179) 

In the off-diagonal one-quark sector {bd and b^d"^) the flow equations are given by 


dSgen(0 

dl 


, r dkdq Dij{k q') 
gJ (27r)® ^ 2uj{k — q) 


(180) 


X Y.^vl{-k)aiVs'{-q)){v\,{-q)ajUs{k)) - Y^{vl{-k)aiUs'{q)){ul,{q)ajUs{k)) 

\ s' / 

X [go{k, q,k- q)gi{k, q,k-q)+ gi{k, q,k- q)go{k, q,k- q)]bs{k)ds{-k ), 


and 

dEj geniP) 

dl 


f dkdq^ D{j(^k q^ 
' 2uj{k - q) 


(181) 
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\ s' s' / 

X [7]o{k, q,k~ q)gi{k, q,k-q)+ qi{k, q,k- q)go{k, q,k- q)]bl{k)dl{-k) , 

where again the symmetry property of go and go, when interchanging any of their two 
arguments, was used. The self-energy operator, at the scale / = 1/A^, is given by 




dkdq ^ 2 Dij{k — q) 


(2i)' 


ft'/ 


2uj{k — q) 


X Y.(yl{-k)aiUs'{q)){u\,{q)ajUs{k)) - Y^{vl{-k)aiVs'{-q)){vl,{-q)ajUs{k)) 




(182) 


and 

= E 


dkdq ^ 2 Dij{k — q) 


\Cf9^ 


(27r)® •' 2uj{k — q) 

X ( Y.iul{k)(^iUs'{q)){ul,{q)ajVs{-k)) - Y^{ul{k)aiVs>{-q)){yl,{-q)ajVs{-k)) 


X 


Do + %-iDl*D^,V>.%Hk)4(_-k) , 


Dg + Df 


(183) 


with the energy denominators dehned in Eq. ( |179|) . Using the spinors, Eq. (|^), and the 
polarization sum, Eq. (|^, the self-energy operators are simplihed to 


So„(A) = E 


dkdq 1 
(27r)6 2 


(1 -h s{k)s{q) + c{k)c{q)k ■ Iq ■ i)^e 


i-o 1 


(1 — s{k)s{q) — c{k)c{q)k ■ Iq ■ i)—e 

Di 


Cfl 

u;{l) 


) 


X {bl{k)bs{k) + dl{-k)d,{-k)) 

T. / (2(-c(fc)s(9) + s(k)c{g)k . iq ■ ()) 


(27r)6 2 

X {bl{k)dl{-k) + d,{-k)bs{k)), 


t 

UJ{1) 

(184) 


with 


Do = -{E{k) + E{q)+u;{l)) 
= E{k) - E{q) - u;{l), 


(185) 


and I = k — q. For large momenta flowing in the loop, one has Dq ~ Di rsj -iE{q) + 
u;{l)) —2uj{q) = —2|qf|. In this limit the self-energy operator, Eq. ( |184| ), is reduced to 


S,„(A) = E 


dkdq 

(^ 


[s{k)s{q) + c{k)c{q)k ■ Iq ■ f) 


^-2D2/a2 
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X {blik)bsik) + dli-k)dsi-k)) 

dkdq 


E 


. ^ (2vr)6 

X {bl{k)dl{-k) + ds{-k)bs{k)), 

where I = k — q, and the energy denominator is 

D = {E{q)+u;{l)). 

This can be simplihed further 


{-c(k)s(q) + s(k)c(q)k ■ Iq ■ () 


( 186 ) 

(187) 






)e 




X (bl(k)b,(k) + dl(-k)d.,(-k)) 

dkdq 


E 


(^—c{k)s{q) + s{k)c{q)k Iq 1 


{E{q)+uj{l))u;{l) 
Cfg^ 

{E{q) + u;{l)Ml) 


)e 


-Aq^lX^ 


s ^ (2vr)6 

X {bl{k)dl{-k) + ds{-k)bs{k)), (188) 

where we rescaled the cut-off, A —>• \/2A, that does not change the result. We regulate 
the self-energy operator arising by normal ordering the instantaneous interactions Eq. 

o 

dkdq 


Einst(A) — ^ 


(s(fc)s(g) -h c{k)c{q)k ■ q) C'/14+c(fc, q)e 


8 - (27r)6 

X ^l{k)bs{k) + dl{-k)ds{-k)] 

dkdq 


E 


{-c{k)s{q) + s{k)c{q)k ■ q) C'jl4+c(fc, q)e 


8 ^ (27r)'^ 

X + ds(-k)6,(fc)], 


(189) 


where the same regulating function is chosen to match the energy denominators. We fol¬ 
lowed the regularization prescription suggested by Zhang and Harindranath [^], where 
the divergent instantaneous terms arising after normal-ordering are regulated using the 
same cut-off function, i.e. an exponent, as in the divergent dynamical terms, while the 
argument in the exponent matches the energy denominator of the interaction. The com¬ 
plete self-energy operator is 


S(A) — Sj„si(A) -|- Sge„(A). 

B.3 Quark condensate 

The second order flow equation, Eq. (p^, for the condensate term is 


(190) 


dOgen{l) _ iV W [ Dij(^k q 


(191) 


dl ' J (27r)® ^ 2ijj{k — q) 

X Y.H{-k)aiUs'{q)){ul,{q)ajVs{-k))2qo{k, q,k- q)go{k, q, fc - q) 1 , 
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where the volume is V = (27r)^5®(0). Integration yields the correction, 60 = 0{l) — 
0 {Iq = 0), where 0{l) is a condensate for the flow parameter /, related to the energy 
scale A by / = 1/A^. The resulting generated condensate term through second order is 


Ogen{X) = N,V 


dkdq ^ 2 Dij{k — q) 


(2vr) 


\Cf9^ 


2u}{k — q) 


X I I , (192) 


which, using the spinors Eq. ( 0 ). is reduced 

Oj,„(A) = -2N^V j ^^(1 + s{k)s{g) + c(k)c(q)k ■ Iq ■ ’ 

(193) 


where I = k — q, and the energy denominator is 


Do = -{E{k) + E{q)+u;{l)). 


(194) 


The regulated condensate terms, generated and instantaneous, can be summarized 


OgenW = -2N,V 


dkdq 

(27r)6 


X 


(1 + s{k)s{q) + c{k)c{q)k ■ Iq ■ i)e 
Cf9^ \ 


-(g+fc+OVA^ 


and, from Eq. (|163|) , 

0,nst{X) = 2N,V 


{E{k) + E{q) + u}{k - q))uj{k - q)) ' 

dkdq 


(195) 


(2ir) 


(1 — s{k)s{q) — c{k)c{q)k ■ q)e CjV{k, q). 


(196) 


The regulating procedure is the same as above (See |^). Note that in the exponential 
factors we have used free dispersion relation, E{q) ~ oj{q) = \q\, valid for large cut-off 
values. The complete radiative correction to the quark condensate, Eq. (154), is 


0(A) — OinstW + OgenW ■ 


(197) 


C Bogoliubov-Valatin transformation and the gap 
equation with double normal ordering 

The Bogoliubov-Valatin transformation for the pairing (BCS) model relates the operators 
b, d which annihilate bare vacuum |0) to a new basis set B, D which annihilate vacuum 
state |fl), containing the quark condensate, 

Bs{k) = ce/2{k)bs{k) - h{s)se/2dl{-k) 

Ds{-k) = ce/2{k)ds{-k) -F h{s)sQ/2bl{k) , (198) 
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where 0 = 0(fc) is the Bogoliubov-Valatin (BCS) angle, and h{s) is the helicity. Similarly 
the transformed quasiparticle spinors are 


Us{k) = ce/2{k)us{k) - h{s)se/2Vs{-k) 

Vs{-k) = ce/2{k)vs{-k) + h{s)se/2Us{k), (199) 

where the bare particle and quasiparticle spinors are dehned, respectively, as 


and 


Us{k) 


Vs{-k) 


1 / + Sa{k)Xs \ 

V2 - Sa{k){(T ■ k)Xs ) ’ 

1 / - s„(fc)(cr • k)T]_s \ 

V2 + Sa{k)T]_s ) 


u,(k) = -l( ). 

V2 - s$(fc)(cr ■ k)xs ) 

v.(-k) = 1 

V2 + s^{k)r]_s ) 


( 200 ) 


( 201 ) 


Here, Xs and rjs are the standard two-component Pauli spinors of a particle and an 
antiparticle, respectively, with t]_s = —i<J 2 Xs- We include masses in the dehnition of 
spinors, therefore the perturbative {a = a{k)) and nonperturbative ($ = *h(A;)) angles 
appear in the sine and cosine as follows 


sin(Q!(A;)) = Sa{k) 


m 


cos(Q!(fc)) 


Ca{k) 


and 


sin(<h(A;)) = s$(fc) 


M{k) 

+ M2(fc) 


, cos(<h(fc)) 


c$(fc) 


k 

\/k‘^ + m'^ ’ 

k 

^fc2 + M2(fc) ’ 


( 202 ) 


( 203 ) 


where m is the bare mass and M{k) is the effective masses. The nonperturbative angle $ 
is related to the Bogoliubov-Valatin angle 0/2, given in Eq. (|198|) , and the perturbative 
angle a, Eq. ( ^02|) , as 


$ = a -f 0 , 


( 204 ) 


So that the following holds 


sq>{k) 

c$(fc) 


v/FT 

k 


m ,, , k 

ce{k) -\- 




y/WT 






-.ce{k) 


m 






fSe{k) 

fSe{k) 
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( 205 ) 





























To obtain the gap equation with double normal ordering we hrst decompose the 
Hamiltonian, Eq. (|^), in the perturbative basis, using (anti)quark b,d operators with 
spinors u,v, Eq. ( p00| ), and the perturbative angle a, Eq. ( |202|) , and then normal order 
H with respect to the perturbative vacuum |0). We obtain the terms summarized in 
Appendix replacing the angle $ by a. Next we perform the Bogoliubov-Valatin trans¬ 
formation from the perturbative to the BCS vacuum, |0) —>■ |f2), expressing the quark 
operators b,d through B,D, Eq. ( p.98|) . Condensate terms do not change, one-body op¬ 
erators transform as 


{bibs + dldg) —>■ cq{BIBs + DIDs) + -|- DgBg) 

{bldl + dsbs) ^ ce{BlDl + DsBs) - h{s)se{BlBs + ^ 1 ^.), ( 206 ) 


and two-body interactions, expressed only by the angle $ using Eq. (|205|) , coincide with 
the interaction terms of Eq. 


). One should be able to combine the angles a and 0 
into $ for the two-body interactions, unless a mistake is done. Now we normal order the 
obtained Hamiltonian with respect to the new vacuum |0) and obtain an additional set of 
terms. In particular, by normal ordering the BV transformed two-body interactions, we 
have the same terms as obtained before Eq. (161) in one-body sector. These terms depend 
only on $. In addition, there are one-body operators obtained by normal-ordering in |0) 
and BV transformed Eq. (|206|) , which depend on a and <h, and also the BV transformed 
kinetic terms, depending only on d). Combining all the terms, we obtain the gap equation 


ks<s>{k) — mc<i,{k) 

= J (c$(fc)[s$(q) - s„(qr)] - s$(fc)[c$(q) - c«(g)]fc ■ q) (207) 

+ J -(^^CfW{k,q) (c$(fc)[s$(q) - s„(qr)] - s$(A;) [c$(qr) - Ca{q)]k ■ Iq ■ /) , 


where the generated terms are also added. One can express this equation in terms of a 
and 0 angles, using Eq. (|205|) , as 


E{k)se{k) = 

+ 

where the self-energy terms 


(S“<^e(fc)ke(9) - 1] + E“ce(fc)se(«) 
E“se(fc)[ce(q) - 1] + E“se{k)se(q)) , (208) 

include the instantaneous and generated contributions 


S — Einst + Sgen • (209) 

The instantaneous self-energies are given only in terms of a 

= {-Ca{k)sa{q) + Sa{k)ca{q)k-q)CfVL+c{k,q) 

= {-Ca{k)co,{q) - s^{k)so,{q)k ■ q)CfVL+c{k, q) 

= {Sa{k)sa{q) -f- Ca{k)ca{q)k ■ q)CfVL+c{k, q) 

= (s„(fc)c„(qr) -c„(fc)s„(qr)fc.q)C'/VL+c(fc,q), (210) 
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and the following changes should be made in these formulas to obtain the generated 
self-energies, S^en; 


and 


k ■ q ^ k iq i, 


( 211 ) 


VL+c(k,q)^W(k,q). 


( 212 ) 


The double normal ordered gap equation, Eqs. (|207|) and (|208|) , is UV finite even 
with the Coulomb potential, Vl+c ~ (A; — alone. For large momenta, q ~ A, 

<h(q) —>• Oi{q) and 0(q) —0, i.e. CQ^q) 1 and SQ^q) —>• 0, and terms coming from 
normal ordering in |r2) and in |0) cancel each other. Indeed, for high UV momenta two 
vacua |0) and Iff) are viewed as the same, giving the same result of normal-ordered 
terms but with opposite signs. However this gap equation is IR inhnite for the conhning 
potential, Vl+c ~ (fc — q)~'^- Cancelation of the leading IR divergence does not happen, 
as it happened with single normal ordering in Iff), Eqs. (|75D and ([761). At fc ~ q the 
r.h.s. of Eq. (|207| ) behaves 


dq 

Wi- 


■CfVn.c(k,q)se(,k) ~ j dq\k 


1-4 


(213) 


and diverges. We conclude that the prescription of double normal ordering leads to 
divergencies for some potentials and is artihcial. 


D Pion mass in the BCS model and in the method 
of flow equations 


In this appendix we check if pion can have zero mass in BCS and flow equation 
approaches. RPA bound state equation, Eq. (|121[) , reads 


(A - B) {k, q)'ip-{q) = Mtlj+{k) 

A{k, q) = 2e(fc) - dqF^^{k, q) 

B{k, q) = -dqF^y{k, q) , (214) 


where tensors F include the instantaneous / and generated G terms, F = I F G.\n. the 
TT channel, using Eq. (jlRjj) for F, 


^xxik) Q) 

ixy{ki q) 

Gxxiki q') 

Gxy{ki q) 


GfVL+c{kq) [ 1 + s{k)s{q) + c{k)c{q)x ] 

GfVi+cikq) [-1 + s{k)s{q) + c{k)c{q)x] 

x{k‘^ -I- — (1 -I- x^)kq 


( 215 ) 


GfW{kq) 

GfW{kq) 


-1 -7 s{k)s{q) + c{k)c{q 
1 -7 s{k)s{q) + c(fc)c(q) 


(fc - q)2 
x{k‘^ + q^) — (1 -I- x‘^)kq 

{k - qf 

















and, from Eq. ( |6^ for e: in the chiral limit m = 0, 

e{k) = kc{k) + J dqCfVL+c{k,q)[s{k)s{q) + c{k)c{q)x] 


j dqCfW{k, q) 

J dqCfVL+c{k,q 


s{k)s{q) + c{k)c{q 

s{q) 


x{k‘^ + — (1 + x^)kq 


s{k) 


+ / dqCfW{k, q) 


{k-q) 

s{q) 

s{k) 


( 216 ) 


where we have used the gap equation with m = 0, Eq. (0. 


k = 


dqCfVn.c c(k)^^-c(q)x 

'rfqW' [c(fc)£M _ - (1 + 


(217) 


L '®(^) ^ qy 

to get In the above formulas the instantaneous interaction, Eq. (^), is given by 

C.V.M = Pi«) 

and the dynamical interactions generated by flow equations, Eq. ( 0 ) and Eq. (|i8|), are 

Cf9^ 


CfW{k,q) = 
CfW{k,q) = 


uj{k - q)[E{q)+uj{k - q)] 
Cf9^ 

E‘^{q) + a;2(fc — q) ’ 


(219) 


where loop momenta in perturbative corrections are large, \k\ |q| (See Appendix 

It is convenient to represent the generated interactions Eq. (|219|) as 


CfW{k,q) = CfU{k,q) -CfSW{k,q) 
-CfW{k,q) = CfU{k,q) -Cf6W{k,q). 


( 220 ) 


with 


CfU{k,q) = 
Cf6W{k,q) = 
Cf6W{k,q) = 


Cf9^ 


uj‘^{k — q) 
Cf9^ 

uj‘^{k — q) 
Cf9^ 


E{q) 


E{q)+uj{k - q) 

E\q) 

E^{q)+uj^{k-. 


( 221 ) 


a;2(fc — q) 

In BCS, with only instantaneous interaction present and hh = hh = 0, we substitute 
Eqs. ( |215|) and (|216| ) into Eq. (p214|) and get 

s{q) 




s(fc) 


ij-{k)-2 J dqCfVL+c{k,q)'ijj_{q) = M'ip+{k). 


( 222 ) 
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As has been found in Ref. choosing the wave function solution as = s{k), the 

l.h.s. of Eq. (|222|) equals zero leading to zero r.h.s. and M = 0. In other words, by this 
choice the kinetic part, ~ 'il^{k), is equal to the potential part (with an opposite sign), 
~ ipiq), thus interaction reduces pion mass to zero. This statement does not depend 
on particular form of the instantaneous potential V, as long as a solution of the gap 
equation exists for V. In particular, for V ~ k/[q^{q^ + 1)], the solution of the gap 
equation is constant, tgik) = s{k)/c{k) ~ / dqs{q)/{q^ + 1) = const, which reduces the 
bound state equation to 


IJ dqCfVL+c{k,q)[^p-{k)-^l)_{q)] = M'ip+{k), 


(223) 


where, in order to have M = 0, we demand ip-ik) = 'ip-{q) for any k, q. This leads to a 
constant wave function, ‘ip-{k) = const, which does not satisfy a normalization condition 
and means a pion bound state is not localized, that is false. 

Moreover, in the original RPA system of equations, Eq. (|118|), the wave functions 


X{k) = -Y{k) = s{k), 


(224) 


give zero l.h.s. and M = 0 in the BCS p (the same type of equation as Eq. (|221D ). 
However, this solution violates the RPA wave function normalization condition, Eq. 

(liioD, 


I dq [X^iq)X{q) - Y*iq)Y{q) ] = j dqr.{q)i’+{q) = 1, 


(225) 


which follows from the commutation relation of the meson creation/annihilation opera¬ 
tors, Eq. ( 101 ). Also, the covariant dispersion law E{P) = \l does not hold in 
this case, which means a breakdown of covariance. As noted by Le Yaouanc et al. |^, 
the model is not covariant since we have adopted an instantaneous interaction. 

Now we calculate the perturbative correction to the BCS solution using flow equa¬ 
tions. We substitute Eqs. (|215| ), ( ^Ri|) and (|22U|) into Eq. ( |214| ), with both instantaneous 


and generated interactions present, and get 

s{q 


2 / dqCfVL+c{k,q) 


s{k) 


'ip-{k)+2 J dqCfU{k,q) 


s{q) 

s{k) 


i^-ik) 


- 2 


dqCfVL+c{k,q)i)_{q)-2 / dqCfU{k, q)'ip_{q) 


(226) 


dq{—Cf6W{k, q)) 


s{k) 


ip-{k) + 2 y dqCfdW{k, q)'tp_{q) = M-^+(fc) . 


Choosing ‘ip-{k) = s{k) the sum of the hrst four terms becomes zero. The rest two terms 
give the pion mass 


M = 2Cfg^ / dqs{q) 


Eiq)[E{q) - uj{k - q)] 


uj{k - q)[E{q)+uj{k - q)][E^{q) + uj^{k - q)] ’ 


(227) 
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at small k {s{k) ~ 1), where E{q) = q in the chiral limit and uj{k — q) = \k — q\. In 
Eq. (p26|) all terms are positive except for the difference, which is negative after angle 
average in the leading order; i.e. 


{E{q) - uj{k - q))^ = {kx - ^(1 - x^) + 0{k^))^ = < 0 , 


(228) 


where x = k ■ q. Therefore the pion mass is negative 


M 



< 0 . 


(229) 


This cannot hold in reality. This false result can be understood because the leading order 
perturbative correction generally lowers a mass of the ground state. Particular feature 
here is that the dynamical kinetic term (self energy), IT, differs from the dynamical 
potential term, W in the two-quark channel, Eq. (p2(]|) . This insures a nonzero (negative) 
correction. Though the flow equation correction has been estimated at small k and it 
is proportional to it cannot be neglected (since the integral diverges at small q). 
Therefore, if the BCS mass is equal to zero which, as discussed before, violates several 
principles, then the flow equations shift this mass to a negative value. This indicates that 
the statement about zero pion mass in the BCS appears to be wrong. We conclude that 
the BCS does not give the zero mass solution for the pion ground state. 

Having shown that the BCS solution Eq. ( 1^24|) violates the normalization condition 
Eq. ( p25|) and is unphysical, Eq. (|229| ), we estimate the vr meson mass when |X| ^ |y| 
in the BCS model. From Eq. (|121|) , {A -|- 5)'?/’+ = 


M = J dkdqi/j^{k){A + B){k,q)i/j+{q) , 


(230) 


where the normalization condition, Eq. (|225|) , has been used. Using Eq. (|214|) for A and 
B and including only the instantaneous part, the mass is given by 


M = 2 dk^pl{k)^p+{k)e{k) - / dkdq'ipX^k) [h^{k, q) + hy{k, q) ] V’+(g) 


(231) 


For the tt meson the instantaneous kinetic term, Eq. (|216|) , can be expressed through the 
instantaneous potential term, Eq. (|215[) , in the chiral limit as 


e{k) = kc{k) + j dq[h^{k,q) + hy{k,q)] . (232) 

Substituting this representation into Eq. ( |231|) the pion mass reads 

M = 2 J dkijl{k)ij+{k)kc{k) (233) 

-h J dkdq^pKk) [I^^{k,q) +I^y{k,q)] [ip+{k) - V’+(q)] 

= 2 J dk'ifj’!^{k)i/j^{k)kc{k) 

+ ‘2 J dkdqxpl{k)CfVL+cik, q) [s{k)s{q) + c{k)c{q)x] [^p+{k) - Mq)] . 
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Here, the first term is always positive and equals 


2 / dkdqipX{k)i)+{k)CfVL+c{k,q) 


- c{k)c{q)x 


> 0 


(234) 


It is the energy of two bare quarks in the transformed BV vacuum. In the second term, the 
kernel averaged over x is positive, since the term dkdq'ip*{k) [I^x + Ixy] V’(fc) represents 
an effective energy of two quarks which is positive (the sum of two quark self energies). 
For any positive operator Q{k,q) > 0, symmetric under interchange of arguments, the 
following holds; 


0 < y dkdq[ip{k) - '0(g)]*Q(fc, q)['i/’(fc) - -0(9)] 

= 2 J dkdq['ljJ{k)*Q{k,q)^jJ{k) - ^jJ{k)*Q{k,q)^jJ{q)] 

= 2 J dkdq'ilj{k)*Q{k, q)['ilj{k) — 'ilj{q)]. (235) 


Therefore the second term in Eq. (|233|) is also positive, provided a positive RPA mass 
of the pion in the BCS model. As shown above, flow equations improve slightly this 
situation and shift the tt meson mass down. However, in the leading order we are unable 
to get exactly zero mass pion even with flow equations. There might be a possibility to 
reach zero mass pion by extending calculations to the higher orders, approching to the 
CO variant result. 


E RPA for the S' and D wave p mesons 

The RPA equations for the p wave function components X(k), Y(k) for L = 0 and L = 2 
states are 


MnX^{k) = 2e{k)X^{k) - 



— 



— 



— 



-MnY^{k) = 2e{k)Y^{k) - 



— 



— 



— 
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MnX^{k) = 2 e{k)X^ 

(k) - 




— 




— 


fytYZ^(k,q)X^(q: 


— 

iyy?r(k,q)Y-(q)- 

fiy-G-,^(k,q)Y^(q) 

-MnY^{k) = 2 e{k)Y^ 

(k) - 




— 


f lZHk.q)X^(l) 


— 




— 

iyy?fik,q)X^iq)- 

fytY-f(k,q)X^(q: 


(236) 


where the instantaneous terms I are 




I°f(k,q) = I, 


jSS 

^xy 


{k,q) 


= lyy {k, q) = CfVi+cik, q)- 


{1 + s{k)){l + s{q)) 


;(1 — s(fc))(l — s{q)){ 4 :x‘^ ~ 1) + ‘ 2 c{k)c{q)x 


lV2, 


I^xxiKq) = I^^^^{Kq)=CfVL^c{k,qy-^{l-s{k)){l-s{q)){x^-l) 
{k, q) = CfVL+c{k, - s{k)){l - s{q)){x^ - 1) 


DS 

yy 

DD 


lyyyq) = iyyk,q) = CfVL+c{yq)- 


-{l + s{k)){l + s{q)){?>x^-l) 


1 


+ ^(1 - s{k)){l - s{q)){x‘^ + 5) + 2 c{k)c{q)x 

6 


= lyx (k, q) = CfVL+c{k, q)- 


il + s{k)){l - s{q)) 


^(1 - s(fc))(l + s{q)) + ‘^c{k)c{q)x 


lyyyq) = iyyyq)=CfVL+c{yq)- 


^(l + s{k)){l - s{qmx^ - 1 ) 


^(1 - s{k)){l + s{q)) - ^c{k)c{q)x 


Ixy{k,q) = Iyy{k,q)=CfVL^c{k,q)- 


3 

2 V 2 


{1 + s{k)){l - s{q)) 


^(1 - s{k)){l + s{q)){ 3 x^ - 1) - ^^c{k)c{q)x 
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-{1 + s{k)){l - s{q)){3x^ - 1) 


+ ^(1 - s(fc))(l + s(qr))(3a;^ - 1) + ‘^c{k)c{q)x 

b 6 


( 237 ) 


and the generated terms G are 




1/1 /111/1 / 11 /i 2(1— 


{k-qf ) 3 

1 


2 / (1 - x^)kq^ 

c{k)c{q) a; + ^ 


(k-q) 


G^-(fc,q) = G^^^^ik,q) = CfW,ik,q)- 


^(l + ,(fc))(l _,(,)) (3 j, 2_ i_ (1 


(fc-qr) 


_ . ,, / (1 — y/2 ^ (1 — x'^)kq 

+ “(1 ~ 'S(^))(l + 'S(g)) [ 2 I + —c{k)c{q) I —2a; + 


6 


(k-q) 


Gi:^^ik,q) = G;;y^ik,q) = GfW,ik,q)- 


V2 

6 


(l + s(fc))(l-s(qr)) 


{k - q )2 
L _ (1 -a; 2 )A ;2 


+ ^(1 - s(fc))(l + s(9)) (^3a;2 “ 1 “ ) + ^c(fc)c(qr) |^-2a; + 

GT{k,q) = G^y^{k,q) 

= GfW,{k,q)l 


(fc-q )2 ^ 
(1 — x‘^)kq^ 
{k-qY ^ 


;(1 + s{k)){l - s{q)) -(1 - 3xY - 77(1 - xY 


k‘^ + 9q‘^ 


+ -(1 - s{k)){l + s{q)) -(1 - 3a;^) - -(1 - a;^) 


6 

+ dk"^ \ 1 


6 


6 


6 


{k - qYj 


+ oC{k)c{q) -Ax - 


{k - qYj 

(1 — x^kq 


(k-q) 


G^4ik,q) = G^y^ik,q)=GfW2ik,q)- 


il + s{k)){l + s{q)) 


+ ^(1 - 5(fc))(l - «(g))(2a;^ - 1) + ^c(fc)c(qr) l^a; ^ 


G^-(A;,q) = G^yYik,q) = GfW2ik,q)- 


{k-qf y 


+ ^(1 - «(*))(! - s(9)) {l + x^- G 2 ) + (21 + 


6 


(k-q) 


GYy^(k,q) = G^Y(k,q) = GfW2(k,q)- 
6 

GYy^(k,q) = G^Y'{k,q) 


6 


(1 + s(fc))(l + s{q)) 2 — 


(fc - qr)2 

(1 — x‘^)kq\ 
(k-q)^ ) 

3(1 — x‘^)k‘^ 


+ 2^(1 - s{k)){l - s{q)) 1^1 + a;2 - c(fc)c(qr) ( 2a; + 


3 


(fc - qr)2 
(1 — xYkq\ 
(k - q)^ ) 


= GfW2{k,q)- 


-(l + .(fc))(l + .(q)) -(9a;^-!)--(!-a;^) 


A ;2 + g 2 

(k-q) 
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.(fc))(l 


s(g)) 



-(1 — x ^) 


P + A 


+ ^c{k)c{q) 



(1 — x^)kq\ 

{k-qy )\ ’ 

(238) 
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